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ABSTRACT 



This thesis developed a data acquisition system to be 
used in conjunction with the 3.5' x 5.0' low speed wind 
tunnel at the Naval Postgraduate School. Interactive 
graphic programs were developed to aid in data acquisiton 
and analysis. In addition, the internal balance that was 
designed by Major Scott Mair and Major Chris Sargent was 
redesigned to correct some problems encountered with the 
drag component. The balance was also instrumented to record 
the pitch and yaw moment components. A calibration rig was 
designed and constructed in order to evaluate the 
interactions of the different components. The equipment 
used and programs developed for data acquisition and 
analysis were adequate. However, balance calibration 
revealed problems with the calibration rig and location of 
the roll component strain gage. Both of these problems will 
have to be corrected before accurate readings can be 
expected from this balance design. 
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I. 



INTRODUCTION 



A . BACKGROUND 

This project used the 3.5' x 5.0' low speed wind tunnel 
at the Naval Postgraduate School to continue the 
aerodynamic study of the effects of helicopter noses and 
tails on drag conducted by Major Mair [Ref. 1] and Major 
Sargent [Ref. 2]. Their work included the construction of 
nine various helicopter configurations, a sting mounted 
support system and a four-degree-of-freedom balance that 
was internal to the model. They also studied the airflow 
around the various configurations with cotton tufting. 

To reduce and analyze the data, Majors Mair and Sargent 
developed several independent computer programs. These 
programs calibrated the balance, corrected the raw data 
from the wind tunnel runs, computed the equivalent flat 
plate area and produced plots of various parameters for 
comparison. However, due to problems encountered with the 
drag component of the balance, their results were 
inconclusive. 

B . GOALS 

The primary goal of this project was to provide the 
students enrolled in the helicopter design class at the 
Naval Postgraduate School with a laboratory type set-up 
that would develop realistic Equivalent Flat Plate Area 
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information for various helicopter configurations. This 
included creating interactive computer programs that the 
students could use with the personal computer located at 
the wind tunnel. 

The secondary goals of this project were to upgrade the 
two-axis internal wind tunnel balance to a six-axis 
internal balance and to develop a calibration rig to 
evaluate the interactions of the balance components. 

Three landing gear configurations were also designed 
and constructed for future analysis of the drag that they 
add to the helicopter. 
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II. APPROACH TO THE PROBLEM 



A. LANDING GEAR DESIGN 

To provide a realistic representation of landing gear 
used with modern-day helicopters, one type of landing gear 
was selected for each nose shape (Figures 2. 1-2. 3). 
Planviews for the three types of landing gear were prepared 
and are included in Appendix A. 

For the attack nose, a skid type of landing gear was 
constructed of aluminum tubing secured to an aluminum 
plate. This type of landing gear is considered a fixed gear 
but was chosen because of its wide use for numerous 
helicopters. For both the smooth nose and blunt nose a 
simulated retractable gear was constructed of wood 
stubwings and model airplane tires. Threaded inserts were 
mounted in the noses and stubwings to allow easy removal of 
the wheel and strut assemblies. This allowed the models to 
be tested in both the clean and dirty configuration. 

B. INTERNAL BALANCE MODIFICATION 

A modified Mair/Sargent balance, Figures 2. 4-2. 7, was 
used for this project. 

To improve the output recorded from the axial 
component, the cuts (see Figure 2.4) alongside that cavity 
were increased by 1/16 of an inch. In addition, the cavity 
itself was squared off thus reducing the curvature of the 
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Figure 2 . 1 



Attack Nose with gear 




Figure 2.2 Smooth Nose with gear 





Figure 2 . 3 Blunt Nose with gear 
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Figure 2.4 Internal Balance with Modifications 
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surfaces in the cavity. These modifications increased the 
flat area upon which the strain gages could be placed. 

An area (see Figure 2.6), similar to that for the 
pitching moment component, was cut to allow recording of 
the yawing moment component. 

For compatability and increased accuracy, the aluminum 
strain gages were removed and replaced by EA-09-062AQ-350 
stainless steel ones. - These gages were smaller and thus 
allowed for better placement within the cavities. To 
provide a longer life, the gages were bonded to the balance 
with an M-Bond AE-15 adhesive system. They were cured at a 
temperature of 150 degrees Fahrenheit for two hours. 

It was desired to record the six component forces on 
the helicopter; lift, drag, yaw, pitching moment, yawing 
moment and rolling moment. However, since the internal 
balance was designed to record only four components, the 
sting support was instrumented to record the yaw force and 
rolling moment (Figure 2.8). 

The gages for the yaw component were placed on the side 
of the sting support to undergo tension and compression 
when subjected to a yawing force. The gages for the rolling 
moment component were placed on top and bottom of the sting 
support at a 45 degree angle to the sting axis. Thus, they 
experienced tension and compression when the model was 
subjected to a rolling moment. 
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Figure 





Figure 2 . 7 Drag Component Gage 




Figure 2.8 Yaw and Roll Moment Gages 



IS 



C. BALANCE CALIBRATION 

In order to determine any interaction between the six 
recording components, a calibration rig was designed and 
mounted to the main fuselage (Figure 2.9). The pans and 
pullies were arranged to simulate the twelve forces and 
moments that the helicopter models would experience in the 
wind tunnel. Each pan, or component was loaded from zero to 
twenty pounds in one pound increments. With each loading, 
readings from the other five components were recorded. This 
method was conducted until all six components were loaded 
in both the positive and negative direction and produced 
thirty interaction matrices. 

A balance calibration program, Figure A. 16, was written 
to determine the relations required to convert the raw data 
counts to actual forces and moments. For each component 
loading, the prime gage constants were determined using the 
following least squares curve fit [Ref. 3]: 

[ Sum ( Xi ) A 2 ] *K1 + [Sum (Xi A 3 ) ] *K2 = Sum(Xi*Yi) 

[Sum(Xi) A 3 ] *K1 + [Sum(Xi A 4) ] *K2 = Sum( (Xi A 2) *Yi) • 

where Xi is the raw data count 
and Yi is the applied load 

Once these constants were computed for all twelve loadings, 
the raw data counts were converted to forces and moments. 
Then for each interaction matrix, the same least squares 
equations were used to determine the interaction 
coefficients. These prime gage constants and interaction 
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Figure 2.9 Balance Calibration Rig 



17 






coefficients were stored for use in the reduction of the 
raw data recorded during the wind tunnel tests. 

D. DATA ACQUISITION 

To improve data reduction, it was desired to have a 
data acquisition system built into the wind tunnel 
instrumentation. The strain gages on' the internal balance 
and support sting were connected to a Pacific Instruments 
model 8255 transducer amplifier (Figure 2.10). The outputs 
from this unit were then routed to data acquisition cards 
mounted in an IBM PC AT. A data acquisition program, Figure 
A. 10, was incorporated into the interactive program set 
that allowed for quick and easy recording of test data. 

Having the transducer amplifier connected to the strain 
gages allowed for easy zeroing and balancing of the gages. 
The data acquisition set-up provided real time readout of 
the forces that the helicopter was experiencing in the wind 
tunnel . 
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Figure 2.10 Data Acquisition Equipment' 
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III. SOLUTION TO THE PROBLEM 



A. DATA COLLECTION AND REDUCTION 

The data collection program, Figure A. 10, was designed 
to record the data in a standardized collection method. 
Interactive steps instruct the user to follow the 
procedures listed below. 



1) . Zero all amplifiers without the model in place. 

2) . Mount the model to the sting support and record 

a zero force reading. 

3) . Place the calibration switches to the + position 

and record a calibration reading. 

4) . Replace the calibration switchres to the center 

position and start the wind tunnel. 

5) . Record first and last data point at zero angle of 

attack with all other readings between plus ten 

degrees and minus eight degrees angle of attack. 

By using a delay loop, changes could be made to the 
angle of attack of the model before the next set of data 
points were recorded. For each tunnel speed, the data 
recorded for the various angles of attack were stored on 
separate files. 

Once the tunnel runs were completed and the data was 
recorded on file, another program, Figure A. 11, was 
designed to convert the raw readings to readings of forces 
and monemts. The sixty interaction coefficients developed 
from the balance calibration were incorporated into this 
data reduction. The following equations, obtained from Mr. 
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David Backs at the NASA AMES Reseasrch Center, were used to 
correct for the interaction of the balance components. 

FI' = FI - (delFl/delF2 ) *F2 - (delFl/delF3 ) *F3 - 

- (delFl/delF2 A 2 ) *F2 A 2 - - (delFl/delF6 A 2 ) *F6 A 2 

Fl"= FI - (delFl/delF2) *F2 ' - 

- (delFl/delF2 A 2) *F2' A 2 - (delFl/delF6 A 2 ) *F6 ' A 2 

These equation were written for each component and 
placed in an iterative loop that checked the difference 
between the two prime values. After the forces were 
corrected for component interaction, a weight tare equation 
was used to correct for the weight of the model 

B. ANALYSIS OF DATA 

The main file of the interactive programs, Figure A. 9 \ 
controlled the data recording, data reduction and data 
analysis with a menu format. For data analysis, the stored 
converted data readings were used to calculate the lift and 
drag coefficients and equivalent flat plate area for each 
angle of attack. Files were created for coefficient of lift 
versus angle of attack, coefficient of drag versus 

coefficient of lift, coefficient of drag versus coefficient 
of lift squared and equivalent flat plate area versus angle 
of attack. There was one file of each created for each 
tunnel speed. 

A plotting routine, Figure A. 15, was included in the 
interactive programs to allow quick analysis of the 
recorded data. To aid in the analysis, up to three plots 
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could be shown on one graph. A delay loop was incorporated 
into the plotting routine to allow for the option of 
obtaining a hard copy of the graph by using the print 
screen command. 
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IV. 



RESULTS 



A. BALANCE CALIBRATION 

The balance was loaded in the twelve component 
directions and the balance calibration program was used to 
produce the tables B.2 through B.13. During the balance 
calibration, a drift in each of the components was noted 
under steady state conditions. To correct for this problem 
it was assumed that each component had a constant drift 
rate. The difference between the first and last zero 
reading was divided by the number of data points taken. 
This correction factor was then applied to each data point. 
This method produced good correlation between different 
data points taken for the same load. 

B. DATA COLLECTION 

Test runs were conducted to evaluate the data 
acquisition program. The vibrations of the model in the 
wind tunnel resulted in erratic fluctuations in the sense 
indicators of the amplifiers. To correct for this problem, 
the acquisition program was modified to collect one hundred 
samples at a rate of five hundred samples per second. The 
average of these one hundred samples was taken as one data 
point. This method produced constant readings for data 
points of similar conditions. 
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C. DATA REDUCTION AND ANALSYIS 

The data reduction program was used to convert the raw 
test data. The interaction equations diverged instead of 
converging to a single value. This pointed out a problem 
with the interaction coefficients. Examination of the 
calibration tables revealed extremely large interactions 
between the loading of the lift component and reaction in 
the roll moment component. There was also noted a large 
interaction between the loading of the pitch moment 
component and reaction in the lift component. 

The reduction program was modified to correct for the 
interactions between the lift and drag components only. The 
reduction program was again executed using the recorded 
test data. This time the interactions converged. This 
proved the validity of the interaction equations and 
confirmed the problems with the balance calibration. 

A test file was created to display the plotting 
capabilities of the interactive programs. Figure 4.1 shows 
the results. 
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TEST FILES FOR PLOTTING 




Figure 4.1 Output from Plotting Routine 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The data acquisition system and programs developed were 
adequate for recording the forces and moments experienced 
in the wind tunnel. The extremely large interactions 
between the lift and roll moment components and between the 
pitch moment and lift components (see Tables B.2 - B.13) 

precluded any analysis of wind tunnel data. However, the 
data acquisition system and balance calibration program 
warrant further research and development of the balance 
system. 

The interactive programs greatly reduced the workload 
required in the data acquisition and analysis phases. They 
provided a quick and easy means for the analysis of the 
recorded data. 

B. RECOMMENDATIONS 

The following are given as recommendations to improve 
the calibration of the internal balance system. 

1 . Balance Modification 

The location of the roll moment strain gage was 
decided to be the cause of the large interaction between 
the lift and roll moment component. This gage was located 
on the back end of the sting support system. The roll 
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moment component was calibrated by applying a torque to the 
center section of the model that was mounted to the 
internal balance. When the lift component was loaded, a 
large bending moment was felt by the roll moment gage and 
thus producing the large interactions. 

By placing the roll moment strain gage at a 
forty-five degree angle on the same cut-out section as the 
pitch moment gage, the torque applied to the center section 
can be used for calibration. Also, the loading of the lift 
component will not greatly affect the reaction of the roll 
moment component. 

2 . Calibration Rig Modification 

By using the center section of the model to mount 
the calibration pans for the pitch moment component, the 
strains were incorrectly transmitted to the lift component. 
A separate calibration set-up will have to be designed that 
isolates the force appplied to the pitch moment component 
from the rest of the balance. 
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APPENDIX A 



SKETCHES AND PROGRAMS 




Figure A.l Nose Gear for Smooth Hose 
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Figure A. 2 Main Gear for Smooth Nose 
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Figure A. 3 Stubwing for Smooth Nose ’(Top View) 
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Figure A. 4 



Skid Gear and Wing for Attack Nose 
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Figure A. 5 Skid Gear for Attack Nose (Top View) 
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Figure A. 6 Nose Gear for Blunt Nose 
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Figure A. 7 Main Gear for Blunt Nose 
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Figure A. 8 Stubwing for Blunt Nose (Top View) 
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REM MAIN. BAS (MAIN CONTROL PROGRAM) 

' PATRICK A. WITT 

' 20 JULY 1985 



REM 

KEY OFF 
SCREEN 1 
COLOR 1 

FOR I = 1 TO 2 

IF I = 1 THEN FILL = 11 ELSE FILL = 12 

IF I = 1 THEN EDGE = 9 ELSE EDGE =10 

IF I = 1 THEN Xl=106 : X2=0: X3=56 ELSE Xl=108: X2=2 : 

X3=58 

IF I = 1 THEN Yl=48: Y2=94: Y3=140 ELSE Yl=46: Y2=92 : 
Y3=13 8 

. .. H i. 

DRAW "C=EDGE ; BM=X1 ? , =Y1 ;U3 0R4D12R12U12R4D3 0L4U12 " 

DRAW "L12D12L4BE1P=FILL; ,=EDGE; " 

i np;ii 

Xl= XI + 26 

DRAW "C=EDGE;BM=X1; , =Y1 ;U3 0R2 0D4L16D9R8D4L8D9R16D4L2 
0BE1P=FILL; , =EDGE ; " 

i ii L n 

Xl= XI + 26 

DRAW "C=EDGE;BM=X1; , =Y1 ;U3 0R4D26R16D4L20BE1P=FILL ; , 
=EDGE ; " 

I IIQII 

Xl= XI + 26 

DRAW "C=EDGE;BM=X1; ,=Y1 ;U30R20D30L20BR4BU4U22R12D22L 
12BG1P=FILL; , =EDGE ; " 

'NEXT LINE 

I IIQII 

DRAW "C=EDGE;BM=X2 ; ,=Y2 ;U30R20D4L16D22R16D4L20BE1P= 
FILL; ,=EDGE; " 

I HQ" 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2 ; ,=Y2 ;U30R20D30L20BR4BU4U22R12D22L 
12 BG1P=FILL ; , =EDGE ; " 

i iipjii 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2 ; ,=Y2 ;U30R4F12U12R4D30L4H12D12L4BE 
1P=FILL; ,=EDGE; " 

i tipii 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2 ; ,=Y2 ; U3 0R2 0D4L16D6R8D4L8D16L4BE1P= 
FILL; ,=EDGE; " 

i ii j ii 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2 ; , =Y2 ;U30R4D30L4BE1P=FILL; , =EDGE ; " 

I II Qll 

Figure A. 9 MAIN. BAS - Main Controlling Program 



36 



640 X2= X2 + 14 

650 DRAW "C=EDGE;BM=X2; ,=Y2;U30R20D4L16D22R12U4L4U4R8" 
660 DRAW "D12L20BE1P=FILL; ,=EDGE; " 

670 ' "U" 

680 X2= X2 + 26 

690 DRAW "C=EDGE ; BM=X2 ; , =Y2 ;U3 0R4D2 6R12U26R4D3 0L2 0BE1P= 
FILL; , =EDGE ; " 

700 ' "R" 

710 X2= X2 + 26 

720 DRAW "C=EDGE ; BM=X2 ; , =Y2 ;U3 0R20D15L12F15L4H15D15L4BU1 
9BR4U7 " 

730 DRAW "R12D7L12BL1P=FILL; ,=EDGE;" 

740 ' "A" 

750 X2= X2 + 30 

760 DRAW "C=EDGE;BM=X2; ,=Y2;U30R20D30L4U15L12D15L4BU19BR 
4U7 " 

770 DRAW "R12D7L12BL1P=FILL; ,=EDGE;" 

780 • "T" 

790 X2= X2 + 26 

800 DRAW "C=EDGE;BH=X2; ,=Y2;BU26U4R20D4L8D26L4U26L8BE1P= 
FILL; ,=EDGE;" 

810 ' "I" 

820 X2= X2 + 26 

830 DRAW "C=EDGE;BM=X2; ,=Y2;U30R4D30L4BE1P=FILL; ,=EDGE;" 
840 ' "0" 

850 X2= X2 + 14 

860 DRAW "C=EDGE;BM=X2; ,=Y2;U30R20D30L20BR4BU4U22R12D22L 
12BG1P=FILL; , =EDGE ; " 

870 ' "N" 

880 X2= X2 + 26 

890 DRAW "C=EDGE;BM=X2; ,=Y2;U30R4F12U12R4D30L4H12D12L4BE 
1P=FILL; , =EDGE ; " 

900 'NEW LINE 
910 ' "A" 

920 DRAW "C=EDGE;BM=X3 ; ,=Y3 ;U30R20D30L4U15L12D15L4BU19BR 
4U7 " 

930 DRAW "R12D7L12BL1P=FILL; , =EDGE ; " 

940 ' "N" 

950 X3= X3 + 26 

960 DRAW "C=EDGE ; BM=X3 ; , =Y3 ;U3 0R4F12U12R4D3 0L4H12D12L4BE 
1P=FILL; , =EDGE ; " 

970 ' ''A'' 

980 X3= X3 + 26 

990 DRAW "C=EDGE;BM=X3 ; ,=Y3 ;U30R20D30L4U15L12D15L4BU19BR 
4U7 " 

1000 DRAW "R12D7L12BL1P=FILL; , =EDGE ; " 

1010 ' "L" 

1020 X3= X3 + 26 

1030 DRAW "C=EDGE;BM=X3;,=Y3;U30R4D26R16D4L20BE1P=FILL;,= 
EDGE ; " 
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1040 ' "Y" 

1050 X3= X3 + 26 

1060 DRAW "C=EDGE;BM=X3 ; ,=Y3 ; BU3 0BL1R5F9R4E9R5G14D16L4 " 
1070 DRAW "U16H14BR2BD1P=FILL; ,=EDGE; " 

1080 ' "S" 

1090 X3= X3 + 40 

1100 DRAW "C=EDGE;BM=X3 ; ,=Y3 ; BU1U2E2R16H19U4E2R2 0F2D2G2 " 
1110 DRAW ''L16F19D4G2L20H2BR2BU1P=FILL; , =EDGE ; " 

1120 ' "I" 

1130 X3= X3 + 32 

1140 DRAW "C=EDGE;BM=X3 ; ,=Y3 ;U30R4D30L4BE1P=FILL; ,=EDGE; " 
1150 ' "S" 

1160 X3= X3 + 14 

1170 DRAW "C=EDGE;BM=X3 ; ,=Y3 ;BU1U2E2R16H19U4E2R20F2D2G2" 
1180 DRAW "L16F19D4G2L20H2BR2BU1P=FILL; ,=EDGE; " 

1190 NEXT I 
1200 LOCATE 22,4 

1210 INPUT "WOULD YOU LIKE INSTRUCTIONS (Y/N)";Z$ 

1220 IF Z$ = "Y» OR Z$ = "y" THEN GOTO 1240 

1230 GOTO 2000 

1240 SCREEN 0 

1250 WIDTH 80 

1260 COLOR 15,1 

1270 CLS 

1280 PRINT 

1290 PRINT TAB ( 15 ) "THIS IS A MENU DRIVEN PROGRAM THAT 
WILL ALLOW YOU TO" 

1300 PRINT TAB (10) "EVALUATE DIFFERENT HELICOPTER CONFIGU 
RATIONS AND DETERMINE" 

1310 PRINT TAB(10) "THEIR LIFT AND DRAG COEFFICIENTS AND 
THEIR EQUIVALENT FLAT" 

1320 PRINT TAB (10) "PLATE AREA. THE DATA GENERATED DURING 
THE TUNNEL RUNS WILL" 

1330 PRINT TAB (10) " BE USED WITH THIS PROGRAM." 

1340 PRINT TAB ( 15) "THE FIRST OPTION IS USED TO RECORD 
THE DATA DURING" 

1350 PRINT TAB (10) "THE TUNNEL RUNS. THE SECOND OPTION 
CONVERTS THE RECORDED" 

1360 PRINT TAB (10) "RAW COUNTS TO FORCES AND MOMENTS. YOU 
WILL HAVE THE CHANCE" 

1370 PRINT TAB ( 10) "TO PROVIDE NAMES FOR EACH OF THESE 
FILES AS YOU GO ALONG." 

1380 PRINT TAB (10) "THE THIRD OPTION LETS YOU EXAMINE THE 
FORCES RECORDED. IT" 

1390 PRINT TAB ( 10 ) "WILL BE OF NO USE TO YOU TO EXAMINE THE 
RAW COUNTS." 

1400 PRINT TAB (15) "THE FOURTH OPTION USES THE CONVERTED 
DATA TO COMPUTE" 

1410 PRINT TAB ( 10 ) "THE DESIRED PARAMETERS. THE DATA 
MUST BE CONVERTED USING" 
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1420 PRINT TAB (10) "OPTION TWO BEFORE THE PARAMETERS CAN BE 
CALCULATED. YOU WILL" 

1430 PRINT TAB (10) "THEN PROVIDE A NAME FOR THE FILES THAT 
STORE THE CALCULATED" 

1440 PRINT TAB (10) "PARAMETERS. USING THESE FILE NAMES, 
THE FIFTH OPTION LETS" 

1450 PRINT TAB (10) "YOU EXAMINE THE DATA THAT CAN BE 
PLOTTED . " 

1455 PRINT: INPUT "PRESS RETURN TO CONTINUE" ;N 

1456 CLS 

1460 PRINT TAB (15) "OPTION SIX MUST BE USED TO SORT THE 
X VALUES BEFORE" 

1470 PRINT TAB ( 10) "THEY CAN BE PLOTTED. THE SEVENTH OPTION 
LETS YOU ADD A GRID" 

1480 PRINT TAB ( 10) "TO THE PLOT AND ALSO MARK THE DATA 
POINTS. THE PLOTTING IS" 

1490 PRINT TAB (10) "IS EXECUTED WITH THE EIGHTH OPTION. 
YOU CAN PLOT UP TO" 

1500 PRINT TAB (10) "THREE CURVES PER PLOT AND RECEIVE A 
HARD COPY BY PRESSING" 

1510 PRINT TAB (10) "THE CTRL-PRTSC BUTTONS. THE FINAL 
OPTION TERMINATES THIS" 

1520 PRINT TAB (10) "PROGRAM. " 

1530 PRINT: INPUT "PRESS RETURN TO CONTINUE" ;N 

2000 ' PLOTTER SET-UP 

2010 SCREEN 0: WIDTH 80: COLOR 14,0,7: CLS: LOCATE 10,1 

2020 PRINT "WHICH DISK DRIVE DO YOU WANT TO STORE YOUR 
DATA FILES ON" ; 

2030 PRINT 

2040 PRINT "DRIVE A IS EITHER ON YOUR LEFT OR ON TOP" 

2050 PRINT "DRIVE B IS EITHER ON YOUR RIGHT OR ON THE 
BOTTOM" 

2070 PRINT 

2080 INPUT "ENTER [ A OR B ]";FD$: FD$=FD$+":" : PD$="C: " 
2090 GRD$="NO GRID": MRK$="MARK " 

2100 0 PT $ = " CURVE " 

2110 GOSUB 2230 
2120 GOSUB 2400 

2130 X%=5+DX%: Y%=12 : LOCATE X%,Y%, 1,0,7 
2140 FIRST%=5: LAST%=13 

2150 ANS$=INKEY$: IF ANS$="" THEN 2150 

2160 IF ANS$=CHR$ (0) +CHR$ (80) THEN IF X%<LAST% THEN X%= 
X%+1 ELSE X%=FIRST% 

2170 IF ANS$=CHR$ ( 0) +CHR$ (72 ) THEN IF X%>FIRST% THEN X%= 
X%-1 ELSE X%=LAST% 

2180 LOCATE X% , Y% : IF ANS$<>CHR$ ( 13 ) THEN 2150 

2190 FLAG1%=CSRLIN: FLAGl%=FLAGl%-4 

2200 ON FLAG1% GOSUB 3120,2740,2630,2800,2680,3170,2870, 
3550,3680 

2210 IF FLAGlIol THEN GOTO 2120 ELSE GOTO 2130 
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222 0 ' INITIALIZING CONSTANTS- 

2230 DIM X(300) ,B(300) , C ( 3 00 ) ,D(300) , Y ( 3 00 ) 
2240 DIM FILESTK$ ( 10) 

2250 OUT 985,6 
2260 SX%=5 : SY%=55 

2270 M$="COMMAND : " 

2280 BLANK$=" 

2290 B$=" " 

2300 SCREEN 0,1: CLS 
2310 HDPOS%=48 

2320 KEY ( 1) ON: ON KEY (1) GOSUB 2350 
2330 RETURN 

2340 ' CHANGE FOREGROUND COLORS 

2350 FG%=FG%+1 

2360 IF FG%>15 THEN FG%=1 

2370 OUT 985, FG% 

2380 RETURN 



2390 


1 


— PRINT MAIN OPTION MENU 




2400 


SCREEN 0,1: 


COLOR 14,0, 


7: CLS: LOCATE 3,1 




2410 


PRINT 


TAB ( 5 ) 


********** 


MAIN OPTIONS 


MENU 




******** *• 








2420 


PRINT 










2430 


PRINT 


TAB (5) 
* fl 


If * 


RECORD TUNNNEL DATA 


2440 


PRINT 


TAB (5) 

* 11 


*** 


REDUCE RAW TUNNEL 


DATA 


2450 


PRINT 


TAB (5) 


*** 


EXAMINE DATA FILE 


FOR TUNNEL 




DATA 










2460 


PRINT 


TAB (5) 
* " 


*** 


CALCULATE Cd, Cl, 


AND E.F.P., 


2470 


PRINT 


TAB ( 5 ) 


II* 


EXAMINE DATA FILE 


FOR PLOTTI 




NG 


*u 








2480 


PRINT 


TAB ( 5 ) 

* ii 


*** 


SORT X VALUES 




2490 


PRINT 


TAB ( 5 ) 

* ii 


II* 


PLOTTING OPTION MENU 


2500 


PRINT 


TAB ( 5 ) 

* *' 


II * 


DO PLOTTING 




2510 


PRINT 


TAB ( 5 ) 

* ii 


II* 


EXIT 




2520 


PRINT 











2530 PRINT TAB (5) "************************************** 

******** ** 



2540 LOCATE 3,65,1: PRINT "STATUS" 

2550 LOCATE 4,55 : PRINT " " 

2560 LOCATE SX%,SY% : PRINT "DATA FILE DRIVE = "+FD$ 
2570 LOCATE SX%+2,SY% : PRINT "USER OPTIONS :" 

2580 LOCATE SX%+3,SY% : PRINT OPT$ 

2590 LOCATE SX%+4,SY% : PRINT GRD$ 

2600 LOCATE SX%+5,SY% : PRINT MRK$ 
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2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 



RETURN 

' DATA EDITOR FOR TUNNEL DATA— 

CHAINFILE$ = PD$+"ADATA. BAS" 

DX%=2 

COMMON FD$ , PD$ , DX% 

CHAIN CHAINFILE$ 

RETURN 

i DATA EDITOR FOR PLOTTING DATA' 

CHAINFILE$ = PD$+ " BDATA . BAS " 

DX%=5 

COMMON FD$ , PD$ , DX% 

CHAIN CHAINFILE$ 

RETURN 

' DATA REDUCTION 

CHAINFILE$ = PD$+"RED. BAS " 

DX%=3 

COMMON PD$ , FD$ , DX% 

CHAIN CHAINFILE$ 

RETURN 

' CALCULATE ROUTINE 

CHAINFILE$ = PD$ + "COMP. BAS" 
DX%=4 

COMMON PD$ , FD$ , DX% 

CHAIN CHAINFILE$ 

REUTRN 

' PLOTTING OPTION MENU 



CLS : 
PRINT 
***'• 


SCREEN 
TAB (20) 


0,1: LOCATE 6,1 

"****** PLOTTING OPTION MENU 


* * * 


PRINT 

PRINT 

*n 


TAB (20) 


If * 


GRID & MARK 




PRINT 

*n 


TAB (2 00 


II * 


GRID & NO MARK 




PRINT 

*m 


TAB (20) 


II* 


NO GRID & MARK 




PRINT 

*» 


TAB (20) 


H * 


NO GRID & NO MARK 




PRINT 

*m 


TAB (20) 


"* 


EXIT 




PRINT 

PRINT 

***" 


TAB (20) 


"************************************* 


X%=8 : 


Y%=27 : 


LOCATE X%, 


Y% , 1,0,7 





FIRST%=8 : LAST%=12 

ANS$=INKEY$ : IF ANS$="" THEN 2990 

IF ANS $=CHR$ ( 0 ) +CHR$ (80) THEN IF X%<LAST% THEN X%= 
X%+1 ELSE X%=FIRST% 

IF ANS$=CHR$ (0) +CHR$ (72) THEN IF X%>FIRST% THEN X%= 
X%+1 ELSE X%=LAST% 
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3020 LOCATE X% , Y% : IF ANS$<>CHR$ ( 13 ) THEN 2990 

3030 FLAG%=CSRLIN: FLAG%=FLAG%-7 

3040 ON FLAG% GOSUB 3070,3080,3090,3100,3110 

3050 DX%=7 

3060 RETURN 

3070 GRD$="GRID MRK$="MARK " : RETURN 

3080 GRD$="GRID " : MRK$="NO MARK" : RETURN 
3090 GRD$="NO GRID": MRK$="MARK " : RETURN 
3100 GRD$="NO GRID": MRK$="NO MARK" : RETURN 
3110 RETURN 

3120 ' RECORD TUNNEL DATA 

3130 CHAINFILE$ = PD$ + "RUNS. BAS" 

3140 DX%=1 

3145 COMMON PD$,FD$,DX% 

3150 CHAIN CHAINFILE$ 

3155 RETURN 

3160 ' SORTING 

3170 CLS : INPUT "ENTER NAME OF FILE TO BE SORTED" ; FILE$ : 
FILE$=FD$+FILE$ 

3180 PRINT "READING FILE "+FILE$+" . . . " : BEEP 
3190 GOSUB 3370 
3200 PRINT "SORTING" 

3210 L%=2 : K%=NOD%-l: R%=NOD% 

3220 WHILE (L%<=R%) 

3230 FOR J%=R% TO L% STEP -1 

3240 IF (X(J%-1) >X(J%) ) THEN SWAP X ( J% ) , X ( J%-1 ) : SWAP 

Y(J%) , Y(J%-1) : K%=J% 

3250 NEXT 

3260 L%=K%~1 

3270 FOR J% = L% TO R% 

3280 IF (X(J%-1) >X(J%) ) THEN SWAP X ( J%) , X ( J%-1) : SWAP 

Y(J%) ,Y(J%-1) : K%=J% 

3290 NEXT 

3300 R%=K%+1 

3310 WEND 
3320 CLS: BEEP 

3330 PRINT "SORTED FILE ";FILE$;" BEING SAVED GOSUB 

3460 

3340 DX%=6 
3350 RETURN 

3360 ' READING A FILE 

3370 OPEN FILE$ FOR INPUT AS #2 

3380 INPUT #2 , NOD% , Yl$ , Xl$ , CONFIG$ , CONF 

3390 INPUT #2 , Q 

3400 FOR J = 1 TO NOD% 



3410 INPUT #2 , Y ( J) , X ( J) 

3420 NEXT 

3430 CLOSE #2 

3440 RETURN 

3450 • SAVE A FILE 
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3460 OPEN FILE$ FOR OUTPUT AS #1 

3470 WRITE #1 , NOD% , Yl$ , Xl$ , CONFIGS , CONF 

3480 WRITE #1,Q 

3490 FOR J = 1 TO NOD% 

3500 WRITE #1,Y(J) ,X(J) 

3510 NEXT 

3520 CLOSE #1 
3530 RETURN 

3540 • INITIALIZING PLOTTER 

3550 CLS: INPUT "NUMBER OF DATA FILE(S) TO BE PLOTTED ON 
THE SAME PLOT =";NOF% 

3560 PRINT "ENTER FILE NAME (S) :":BEEP 
3570 FOR 1% = 1 TO NOF% 

3580 PRINT "FILE #"? -.PRINT USING "##";I%; : INPUT " =" ? 

FILES 

3590 FILES = FD$ + FILES 

3600 FILESTK$ ( 1% ) = FILES 

3601 PRINT "HOW DO YOU WANT THIS CURVE PLOTTED ? ENTER 
(POINT) TO JUST PLOT" 

3602 INPUT "THE POINTS OR (CURVE) TO PRODUCE A CURVE 
FIT" ;OPT$ 

3610 CURVES (1%) = OPTS 
3620 NEXT 

3630 PFILES = PD$+"MAIN. BAS" : FILE2$= PD$ + "PLOTTER. BAS" 

3640 SL% = 100 

3650 CHAIN FILE2 $ , 3 0 , ALL 

3660 DX%=8 

3670 RETURN 

3680 CLS: 

3690 LOCATE 10,5: PRINT "PROGRAM TERMINATED, REMOVE DATA 
FILES FROM DISK DRIVE" 

3700 LOCATE 11,5: PRINT "AND SECURE COPMUTER AND AMPLIFIE 
S'" 

3710 END 



Figure A. 9 MAIN. BAS (cont.) 



43 



10 
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460 
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480 



NAME: Data Acquisition And Control ( DAAC) 

HEADER for BASICA 

FILE NAME: DACHDR . BAS 

DOS DEVICE NAME: DAAC 

RESERVED FUNCTION NAMES: 

AINM, AINS, AINSC, AOUM, AOUS , 

BINM, BINS, BITINS, BITOUS , BOUM, BOUS , 
CINM, CINS, CSET , DELAY 
RESERVED DEF SEG VALUE NAME: DSEG 

NAMES DEFINED AND USED BY HEADER: 

ADAPT % , AI, COUNT, FOUND% , 

HNAME$ , SG% , STAT% 



When using the BASICA Interpreter, this header 
must be executed before any function calls are 
made that access the DAAC adapter. It initializes 
a number of variables for each function call. These 
variables are reserved and should not be used except 
to access the DAAC adapter. This routine also does a 
DEF SEG to the segment where the DAAC Device Driver 
(DAC.COM) is loaded. If you execute a DEF SEG to 
access other hardware, you must DEF SEG to the segment 
of the DAAC Device Driver before any subsequent 
calls to access the DAAC adapter. 



FOUND% = 0 
SG% = &H2E 

'Start searching the interrupt vectors until you find 
'one that points to the DAAC device driver. 

'Do a DEF SEG to that segment. 

WHILE ( (SG% <= &H3E) AND (FOUND% = 0)) 

DEF SEG = 0 

DSEG = PEEK ( SG% ) + PEEK(SG% + 1) * 256 

DEF SEG = DSEG 
HNAME$=" " 

FOR AI=10 TO 17 

HNAME$ = HNAME$ + CHR$ ( PEEK (AI ) ) 

NEXT AI 

IF HNAME$ = "DAAC " AND PEEK (18) + PEEK (19) <> 
0 THEN FOUND% = 1 
SG% = SG% + 4 



WEND 

IF FOUND% = 
NOT FOUND" : 



0 THEN PRINT "ERROR: DEVICE DRIVER DAC.COM 
END 
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490 'Now initialize all function name variables for calls 
500 'to access the device driver. 



510 


AINM 


= 


PEEK ( &H13 ) 


* 


256 


+ 


PEEK ( &H12 ) 


520 


AINS 


= 


PEEK(&H15) 


* 


256 


+ 


PEEK ( &H14 ) 


530 


AINSC 


= 


PEEK ( &H17 ) 


* 


256 


+ 


PEEK ( &H16) 


540 


AOUM 


= 


PEEK ( &H19 ) 


* 


256 


+ 


PEEK ( &H18 ) 


550 


AOUS 


= 


PEEK(SHIB) 


* 


256 


+ 


PEEK ( &H1A) 


560 


BINM 


= 


PEEK(ScHlD) 


* 


256 


+ 


PEEK ( &H1C) 


570 


BINS 


= 


PEEK(SHIF) 


* 


256 


+ 


PEEK ( &H1E) 


580 


BITINS 


= 


PEEK ( &H21) 


* 


256 


+ 


PEEK ( &H2 0 ) 


590 


BITOUS 


= 


PEEK ( &H23 ) 


* 


256 


+ 


PEEK ( &H22 ) 


600 


BOUM 


= 


PEEK( &H25) 


* 


256 


+ 


PEEK(&H24) 


610 


BOUS 


= 


PEEK( &H27) 


* 


256 


+ 


PEEK ( &H2 6 ) 


620 


CINM 


* = 


PEEK ( &H29 ) 


* 


256 


+ 


PEEK ( &H2 8 ) 


630 


CINS 


= 


PEEK ( &H2B) 


* 


256 


+ 


PEEK ( &H2A) 


640 


CSET 


= 


PEEK ( &H2D) 


* 


256 


+ 


PEEK ( &H2C) 


650 


DELAY 


= 


PEEK ( &H2F) 


* 


256 


+ 


PEEK ( &H2E) 


660 


' Finally, 


execute any call 


to : 


re- 


-initialize the 



670 'device driver from any former invocation of BASIC. 
680 ADAPT% = 0 
690 COUNT = 1 • 

700 STAT% = 0 

710 CALL DELAY ( ADAPT % , COUNT, STAT%) 

720 ' 

730 'End of DAAC BASICA Header 
740 ' 

750 REM — RUNS. BAS; PROGRAM TO RECORD THE DATA FROM 

WIND 

7 60 REM— TUNNEL RUNS 

770 KEY OFF: COLOR 15,1,4: CLS 

780 CLS 

790 PRINT: PRINT 

800 PRINT TAB (10) "HELO CONFIGURATIONS" 



810 


PRINT : 


PRINT 








820 


PRINT 


TAB (5) 


" 1 . 


ATTACK NOSE, 


STRAIGHT TAIL" 


830 


PRINT 


TAB (5) 


" 2 . 


ATTACK NOSE, 


LOW TAIL" 


840 


PRINT 


TAB ( 5 ) 


"3. 


ATTACK NOSE, 


HIGH TAIL" 


850 


PRINT 


TAB ( 5 ) 


"4. 


SMOOTH NOSE, 


STRAIGHT TAIL" 


860 


PRINT 


TAB (5) 


"5. 


SMOOTH NOSE, 


LOW TAIL" 


870 


PRINT 


TAB ( 5 ) 


" 6 . 


SMOOTH NOSE, 


HIGH TAIL" 


880 


PRINT 


TAB ( 5 ) 


"7 . 


BLUNT NOSE, 


STRAIGHT TAIL" 


890 


PRINT 


TAB ( 5 ) 


" 8 . 


BLUNT NOSE, 


LOW TAIL" 


900 

910 


PRINT 
PRINT : 


TAB ( 5 ) 

: PRINT 


"9 . 


BLUNT NOSE, 


HIGH TAIL" 


920 

930 


INPUT 

PRINT 


"WHICH 


CONFIGURATION IS 


BEING RUN" ;N 



940 INPUT "DOES CONFIGURATION INCLUDE LANDING GEAR";Z$ 
950 ON N GOSUB 4550,4620,4690,4750,4820,4870,4920,4990, 
5040 

955 i WEIGHTS ARE IN POUNDS 
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960 


IF 


CONF =10 


THEN 


W = 


22.7892 


970 


IF 


CONF = 15 


THEN 


W = 


23 . 4623 


980 


IF 


CONF =20 


THEN 


W = 


24.6945 


990 


IF 


CONF =25 


THEN 


W = 


25.3676 



1000 


IF 


CONF 


= 


30 


THEN 


W 


= 


24.6945 


1010 


IF 


CONF 


= 


35 


THEN 


W 


= 


25.3676 


1020 


IF 


CONF 


= 


40 


THEN 


W 


= 


20.8134 


1030 


IF 


CONF 


= 


45 


THEN 


W 


= 


21.5866 


1040 


IF 


CONF 


= 


50 


THEN 


W 


= 


22 . 7187 


1050 


IF 


CONF 


= 


55 


THEN 


W 


= 


23.4913 


1060 


IF 


CONF 


= 


60 


THEN 


W 


= 


22 . 7187 


1070 


IF 


CONF 


= 


65 


THEN 


W 


= 


23.4913 


1080 


IF 


CONF 


= 


70 


THEN 


W 


= 


21.3315 


1090 


' IF 


CONF 


= 


75 


THEN 


W 


= 


23 .2892 


1100 


IF 


CONF 


= 


80 


THEN 


W 


= 


23.2368 


1110 


IF 


CONF 


= 


85 


THEN 


W 


= 


25.1945 


1120 


IF 


CONF 


= 


90 


THEN 


W 


= 


23 .2368 


1130 


IF 


CONF 


= 


95 


THEN 


W 


= 


25.1945 



1140 CLS 

1150 REM RECORD TUNNEL DATA — 

1160 DIM L(100) ,D(100) , Y(100) ,PM(100) ,YM(100) ,RM(100) ,AOA( 
100 ) 

1170 DIM DAT(399) , DAT% (399) , DAT1 (399) ,DAT1% (399) 

1180 COLOR 15,1: KEY OFF: CLS 

1190 PRINT "RECORDING OF WIND TUNNEL RAW DATA" 

1200 INPUT "WHAT IS THE TUNNEL SPEED (Q) FOR THIS RUN";Q' 
1210 PRINT 

2670 PRINT "WITHOUT THE MODEL IN PLACE, ADJUST THE ZERO 
SET SCREWS TO ZERO" 

2671 PRINT "OUT EACH AMPLIFIER. AFTER ALL AMPLIFIERS ARE 
ZEROED, LOAD THE" 

2672 PRINT "MODEL ONTO THE STING. WHEN THE MODEL IS 
MOUNTED, PRESS RETURN" 

2673 INPUT "TO RECORD A NO FORCE ZERO READING" ;N 

2690 PRINT: PRINT " ZD ZL ZY ZPM ZYM 

ZRM ZAOA" 

2700 STAT%=0: MODE%=0 : STOR%=0: COUNT=100: RATE=500 

2710 ADAPT %= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 

2720 CALL AINSC ( ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

2730 ZD = 0:ZPM=0:ZL=0:ZYM=0 

2740 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

2750 FOR J = 0 TO 396 STEP 4 
2760 DAT ( J) = ( DAT% (J)/204.8)-10 
2770 ZD =ZD + DAT ( J ) 

2780 NEXT J 

2790 ZD = ZD/100 

2800 FOR J = 1 TO 397 STEP 4 

2810 DAT ( J) = ( DAT% (J)/204.8)-10 
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2820 ZL = ZL + DAT ( J) 

2830 NEXT J 

2840 FOR J = 2 TO 398 STEP 4 
2850 DAT ( J) = (DAT% (J)/204.8)-10 
2860 ZPM = ZPM + DAT ( J ) 

2870 NEXT J 

2880 FOR J = 3 TO 399 STEP 4 
2890 DAT ( J) = (DAT% (J)/204.8)-10 
2900 ZYM = ZYM + DAT ( J ) 

2910 NEXT J 

2920 STAT%=0: MODE%=0 : STOR%=0 : COUNT=100: RATE=500 

2930 ADAPT %= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=2 

2940 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% ,MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

2950 ZY =0 : ZRM=0 : ZAOA=0 

2960 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

2970 FOR J - 0 TO 297 STEP 3 
2980 DAT1 ( J) =(DAT1% (J)/204.8)-10 
2990 ZRM = ZRM + DAT1 ( J) 

3000 NEXT J 

3010 ZRM = ZRM/100 

3020 FOR J = 1 TO 298 STEP 3 

3030 DAT1 (J) =(DAT1% (J)/204.8)-10 

3040 ZY = ZY + DAT1(J) 

3050 NEXT J 

3051 FOR J = 2 TO 299 STEP 3 

3052 DAT1 ( J) = (DAT1% (J)/204.8)-10 

3053 ZAOA = ZAOA + DAT1(J) 

3054 NEXT J 

3060 ZL=ZL/100 : ZPM=ZPM/100 : ZYM=ZYM/100 : ZY=ZY/100 : ZAOA=ZAOA 
/100 

3070 LOCATE 11,1: PRINT USING "+#.###";ZD: LOCATE 11,10: 
PRINT USING "+#.###" ;ZL 
3080 LOCATE 11,19: PRINT USING "+#.###" ;ZY 
3090 LOCATE 11,28: PRINT USING "+#.###" ;ZPM 
3100 LOCATE 11,37: PRINT USING "+#.###";ZYM 

3110 LOCATE 11,46: PRINT USING "+#.###"; ZRM 

3111 LOCATE 11,55: PRINT USING "+#.###"; ZAOA 
3120 REM 

3130 PRINT: PRINT 

3140 INPUT "AFTER PLACING ALL CAL SWITCHES TO + SETTING HIT 
RETURN" ;X 
3150 PRINT 

3160 PRINT " CALD CLL CALY CALPM CALYM 

CALRM CALAOA" 

3170 STAT%=0 : MODE%=0: STOR%=0 : COUNT=100: RATE=500 

3180 ADAPT %= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 

3190 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 
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3200 CALD=0 : CLL=0 : CALYM=0 : CALPM=0 

3210 IF STAT%<> 0 , THEN PRINT USING "EXECUTION ERROR ###" 
;STAT% : END 

3220 FOR J = 0 TO 396 STEP 4 
3230 DAT ( J) = ( DAT% (J)/204.8)-10 
3240 CALD =CALD + DAT ( J ) 

3250 NEXT J 

3260 CALD = CALD/100 

3270 FOR J = 1 TO 397 STEP 4 

3280 DAT ( J ) = ( DAT% (J)/204.8)-10 

3290 CLL = CLL + DAT ( J) 

3300 NEXT J 

3310 FOR J = 2 TO 398 STEP 4 
3320 DAT ( J) = ( DAT% (J) /204 .8) -10 
3330 CALPM = CALPM + DAT ( J ) 

3340 NEXT J 

3350 FOR J = 3 TO 399 STEP 4 
3360 DAT ( J) = ( DAT% (J)/204.8)-lG 
3370 CALYM = CALYM + DAT ( J) 

3380 NEXT J 

3390 STAT%=0: MODE%=0: STOR%=0 : COUNT=100: RATE=500 

3400 ADAPT%= 1; DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=2 

3410 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% ,MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

3420 CALY =0 : CALRM=0 : CALAOA=0 

3430 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

3440 FOR J = 0 TO 297 STEP '3 
3450 DAT1 ( J) = (DAT1% (J)/204.8)-10 
3460 CALRM = CALRM + DAT1(J) 

3470 NEXT J 

3480 CALRM = CALRM/100 
3490 FOR J = 1 TO 298 STEP 3 
3500 DAT1(J)=(DAT1% (JJ/204.8) -10 
3510 CALY = CALY + DAT1(J) 

3520 NEXT J 

3521 FOR J = 2 TO 299 STEP 3 

3522 DATl(J)=(DATl%(J)/204.8) -10 

3523 CALAOA = CALAOA + DAT1(J) 

3524 NEXT J 

3530 CLL=CLL/100 : CALPM=CALPM/100 : CALYM=C ALYM/ 100 : CALY=CALY 
/100 : CALAOA=CALAOA/100 
3540 LOCATE 18,1: PRINT USING "+#.###"; CALD 
3550 LOCATE 18,10: PRINT USING »+#.###";CLL 
3560 LOCATE 18,19: PRINT USING "+#.###"; CALY 
3570 LOCATE 18,28: PRINT USING "+#.###"; CALPM 
3580 LOCATE 18,37: PRINT USING "+#.###"; CALYM 

3590 LOCATE 18,46: PRINT USING "+#.###"; CALRM 

3591 LOCATE 18,55: PRINT USING "+#.###"; CALAOA 

3592 INPUT "PRESS RETURN TO CONTINUE" ;N 

Figure A. 10 RUNS. BAS (cont.) 
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CLS 

PRINT: PRINT "REPLACE THE CAL SWITCHES TO THE CENTER 
POSITION" 

PRINT "AFTER YOU HAVE GOTTEN THE WIND TUNNEL UP TO 
SPEED AND ARE READY" 

PRINT "TO RECORD DATA PRESS RETURN. ONCE THE HEADINGS 
ARE PRINTED THE" 

PRINT "F2 KEY WILL RECORD THE DATA. THE FI KEY WILL 
SAVE THE DATA AND " 

PRINT "RERUN YOU TO THE MAIN MENU. YOU CAN TAKE READ 
INGS FOR ANY ANGLE" 

PRINT "OF ATTACK BETWEEN +10 AND -8 DEGREES." 

PRINT 

PRINT "THIS PROGRAM CORRECTS FOR DRIFT IN THE BALANCE 
AND AMPLIFIERS." 

PRINT "YOUR FIRST AND LAST DATA POINTS SHOULD BE TAKEN 
AT ZERO ANGLE OF" 

PRINT "ATTACK. ALSO, DO NOT TAKE MORE THAN ONE DATA 
POINT FOR THE SAME " 

PRINT "ANGLE OF ATTACK EXCEPT FOR YOUR FIRST AND LAST 
DATA POINT. YOU" 

PRINT "ARE LIMITED TO 100 -DATA POINTS. REMEMBER, EACH 
RUN IS FOR ONE" 

PRINT "SPECIFIC ' Q • SETTING" 

INPUT "PRESS RETURN TO CONTINUE" ;N 
CLS 

PRINT " DRAG LIFT YAW PITCH YAW 



AOA 



LIFT YAW 

TUNNEL" 



MOM. 



MOM. 



SPEED" 



'SET STOP FLAG 
'RECORD DATA 



ROLL 
PRINT " 

MOM. 

SOAP = 0: N=4 
FOR K = 1 TO 100 
ON KEY ( 1 ) GOSUB 3820 
ON KEY (2) GOSUB 3850 
KEY ( 1 ) ON: KEY ( 2 ) ON 
IF SOAP = 2 THEN GOTO 3780 
IF SOAP = 1 THEN GOTO 4370 
GOTO 3720 
SOAP = 0 
NOD% = K 
NEXT K 
GOTO 4370 

REM SET STOP FLAG 

SOAP = 1 

RETURN 

REM STEPS TO RECORD DATA 

STAT%=0: MODE%=0 : STOR%=0: COUNT=100: RATE=500 

ADAPT%= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 
CALL AINSC ( ADAPT % , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 
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3890 D(K)=0:L(K)=0: YM(K) =0 : PM(K) =0 

3900 IF STAT%<> 0 , THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

3910 FOR J = 0 TO 396 STEP 4 
3920 DAT ( J) = (DAT% (J)/204.8)-10 
3930 D (K) =D(K) + DAT ( J) 

3940 NEXT J 

3950 D (K) = D(K)/100 

3960 FOR J = 1 TO 397 STEP 4 

3970 DAT ( J) = ( DAT% (J)/204.8)-10 

3980 L(K) = L(K) + DAT ( J) 

3990 NEXT J 

4000 FOR J = 2 TO 398 STEP 4 
4010 DAT ( J ) = ( DAT% (J)/204.8) -10 
4020 PM (K) = PM (K) + DAT ( J ) 

4030 NEXT J 

4040 FOR J = 3 TO 399 STEP 4 

4041 DAT ( J) = (DAT% (J)/204.8) -10 
4050 YM (K) = YM (K) + DAT ( J ) 

4060 NEXT J 

4070 STAT%=0: MODE%=0 : STOR%=0: COUNT=100: RATE=500 

4080 ADAPT %= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=2 

4090 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% ,MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

4100 Y (K) =0:RM(K)=0: AOA(K)=0 

4110 ‘IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

4120 FOR J = 0 TO 297 STEP 3 
4130 DAT1 ( J) = (DAT1% (J)/204.8)-10 
4140 RM (K) = RM (K) + DAT1(J) 

4150 NEXT J 

4160 RM (K) = RM (K) /100 
4170 FOR J = 1 TO 298 STEP 3 
4180 DAT1 ( J) = (DAT1% (J)/204.8)-10 
4190 Y (K) = Y(K) + DAT1(J) 

4200 NEXT J 

4210 FOR J = 2 TO 299 STEP 3 
4220 DAT1 ( J) = ( DAT1% (J)/204.8)-10 
4230 AOA(K) = AOA(K) + DAT1(J) 

4240 NEXT J 

4250 AOA(K)=AOA(K)/100:L(K)=L(K)/100:PM(K)=PM(K)/100:YM(K) 
=YM(K)/100: Y(K) =Y(K)/100 

4255 B = AOA(K) - ZAOA 

4256 AOA(K) = (5.8469*B) + (. 0077583* (B A 2) ) 

4260 IF N> 23 THEN N=4 : CLS : PRINT " DRAG LIFT YAW 

PITCH YAW ROLL AO A TUNNEL" : PRINT " 

MOM. MOM. MOM. SPEED" 

4270 LOCATE N,l: PRINT USING "+#.###"; D (K) 

4280 LOCATE N,10: PRINT USING "+#.###"; L (K) 

4290 LOCATE N,19: PRINT USING "+#.###"; Y (K) 

Figure A. 10 RUNS. BAS (cont.) 
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4300 LOCATE N,28: PRINT USING "+#.###"; PM (K) 

4310 LOCATE N,37: PRINT USING "+#.###" ?YM(K) 

4320 LOCATE N,46: PRINT USING "+#.###" ;RM(K) 

4330 LOCATE N,55: PRINT USING " + #.###'• ;AOA(K) 

4340 N=N+1 
4350 SOAP = 2 
4360 RETURN 

4370 SOAP = 0 

4371 ' CORRECT FOR DRIFT 

4372 DIFFL =(L(NOD%) -L(l) ) : DIFFD= (D (NOD%) -D ( 1) ) : DIFFY= (Y ( 
NOD%) -Y(l) ) 

4373 DIFFPM=(PM(NOD%) -PM(1) ) : DIFFYM= ( YM (NOD%) -YM ( 1) ) : DIFFR 
M=(RM(NOD%)-RM(l) ) : DIFFAOA= (AOA(NOD%) -AOA(l) ) 

4374 D=N0D%-1 

4375 CORL=DIFFL/D : CORD=DIFFD/D : CORY=DIFFY/D : CORPM=DIFFPM/D 
: CORYM=DIFFYM/D : CORRM=DIFFRM/D : CORAOA=DIFFAOA/D 

4376 A=1 

4377 FOR K = 2 TO NOD% 

4378 L (K) =L (K) - (A*CORL) 

4379 D (K) =D (K) - (A*CORD) 

4380 Y(K)=Y(K) -(A*CORY) 

4381 PM ( K) =PM ( K) - ( A*CORPM) 

4382 YM (K) =YM (K) - (A*CORYM) 

4383 RM (K) =RM (K) - (A*CORRM) 

4385 A = A+l 

4386 NEXT K 

4387 ' SAVE DATA TO FILE 

4390 CLS : PRINT 

4400 INPUT "WHAT IS THE NAME FOR THE FILE TO STORE THE RAW 
DATA" ;FILE$ 

4410 FILE$=FD$+FILE$ 

4420 OPEN FILE$ FOR OUTPUT AS #1 

4421 WRITE # 1 , NOD% , CONFIG$ ; CONF 
4430 WRITE #1,Q,W 

4470 WRITE #1 , ZL, ZD , ZY , ZPM, ZYM, ZRM, ZAOA 

4480 WRITE # 1 , CLL , CALD , CALY , CALPM , CALYM , CALRM , CALAOA 

4490 FOR J=1 TO NOD% 

4500 WRITE #1,L(J) ,D(J) ,Y(J) , PM ( J ) , YM ( J ) , RM ( J ) , AO A ( J ) 

4510 NEXT J 

4520 CLOSE #1 

4521 CLS: PRINT "DO YOU WANT TO TAKE ANOTHER SET OF DATA 
POINTS FOR THE SAME" 

4522 INPUT "CONFIGURATION BUT A DIFFERENT ' Q ' SETTING (Y/N 
) " ; ANS$ 

4523 IF ANS$="Y" OR ANS$="y" THEN GOTO 4524 ELSE GOTO 4530 

4524 PRINT: INPUT "WHAT IS THE ' Q ' VALUE FOR THIS RUN";Q 

4525 GOTO 3610 

4530 ' RETURN TO MAIN PROGRAM 

4531 COMMON PD$,FD$,DX% 

4532 CHAIN PD$+"MAIN . BAS" , 2 110 

Figure A. 10 RUNS. BAS (cont.) 
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REM SUBROUTINES FOR CONFIGURATION MARKING- 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4590 
CONF =10 

CONFIG$= "ATTACK NOSE/STRAIGHT TAIL WITHOUT GEAR" 

GOTO 4610 
CONF = 15 

CINFIG$=" ATTACK NOSE/STRAIGHT TAIL WITH GEAR" 

RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4660 
CONF =20 

CONFIG$=" ATTACK NOSE/LOW TAIL WITHOUT GEAR" 

GOTO 4680 
CONF =25 

CINFIG$=" ATTACK NOSE/LOW TAIL WITH GEAR" 

RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4730 
CONF =30 

CONFIG$=" ATTACK NOSE/HIGH TAIL WITHOUT GEAR" 

GOTO 4740 

CONF = 35: CINFIG$=" ATTACK NOSE/HIGH TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4790 
CONF =40 

CONFIG$=" SMOOTH NOSE/STRAIGHT TAIL WITHOUT GEAR" 

GOTO 4810 
CONF =45 

CINFIG$=" SMOOTH NOSE/STRAIGHT TAIL WITH GEAR" 

RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4850 
CONF=50:CONFIG$=" SMOOTH NOSE/LOW TAIL WITHOUT GEAR" 
GOTO 4860 

CONF = 55: CINFIG$=" SMOOTH NOSE/LOW TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4900 
CONF=60:CONFIG$=" SMOOTH NOSE/HIGH TAIL WITHOUT GEAR" 
GOTO 4910 

CONF = 65: CINFIG$=" SMOOTH NOSE/HIGH TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4960 
CONF =70 

CONFIG$=" ATTACK NOSE/STRAIGHT TAIL WITHOUT GEAR" 

GOTO 4980 
CONF =75 

CINFIG$=" ATTACK NOSE/STRAIGHT TAIL WITH GEAR" 

RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 5020 

CONF = 80 : CONFIG$=" BLUNT NOSE/LOW TAIL WITHOUT GEAR" 

GOTO 5030 

CONF = 85: CINFIG$=" BLUNT NOSE/LOW TAIL WITH GEAR" 
RETURN 
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IF Z$ = "Y" OR Z$ = "y" THEN GOTO 5070 

CONF=90 : CONFIG$=" BLUNT NOSE/HIGH TAIL WITHOUT GEAR" 

GOTO 5080 

CONF = 95: CINFIG$=" SMOOT NOSE/HIGH TAIL WITH GEAR" 
RETURN 
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(CONVERT RAW COUNTS TO FORCES)- 



REM RED . BAS 

REM 
REM 

COLOR 15,1: KEY OFF: CLS 

DIM L ( 100) , LI (100) ,L2 (100) ,D(100) ,D1(100) ,D2 (100) , 

(100) , Y1 ( 100) ,Y2 (100) 

DIM PM (100) , PM1 (100) ,PM2 (100) ,YM(100) ,YM1(100) , YM2 (10 
0) ,RM(100) , RM1 ( 100 ) ,RM2 (100) 

DIM LF(100) , DF(100) , AOA(IOO) 

PRINT : PRINT "WHAT IS THE NAME OF THE FILE THAT CONTAINS 
THE RAW DATA" 

INPUT "THAT YOU WISH TO CONVERT" ; FILE $ 

FILE$=FD$+FILE$ 

' READ IN RAW DATA FILE — 

OPEN FILES FOR INPUT AS #1 
INPUT #1 , NOD% , CONFIG$ , CONF 
INPUT #1, Q, W 

INPUT #1, ZL, ZD, ZY, ZPM, Z YM , ZRM , ZAOA 
INPUT # 1 , CLL , CALD , CALY , CALPM , CALYM , CALRM , CALAOA 
FOR J = 1 TO NOD% 

INPUT #1 , L ( J) , D ( J) , Y ( J) , PM ( J ) , YM ( J ) ,RM(J) ,AOA(J) 
NEXT J 
CLOSE #1 

PRINT: PRINT "PLEASE WAIT WHILE PERFORMING CONVERSION" 
PRINT "THE CONVERSION WILL TAKE APPROXIMATELY TEN MIN 
UTES FOR " 

PRINT "TWENTY DATA POINTS" 

i READ IN CALIBRATION CONSTANTS 

OPEN "C: CONST" FOR INPUT AS #1 

INPUT # 1 , INCALL , INCALD , INCALY , INCALPM , INCALYM , INCALRM 
INPUT #1 , K1LPOS , K2LPOS , K1DPOS , K2DPOS , K1YPOS , K2YPOS 
INPUT # 1 , K1PMPOS , K2 PMPOS , K1 YMPOS , K2 YMPOS , K1RMPOS , 
K2RMPOS 

INPUT # 1 , K1LNEG , K2LNEG , K1DNEG, K2DNEG , K1YNEG , K2YNEG 
INPUT # 1 , K1PMNEG , K2PMNEG , K1YMNEG , K2 YMNEG , K1RMNEG , 
K2RMNEG 

INPUT # 1 , DDDL1P , DDDL2P, DYDL1P , DYDL2 P , DPMDL1P , DPMDL2 P , 

DYMDL1P, DYMDL2P 

INPUT # 1 , DRMDL1 P , DRMDL2 P 

INPUT #1, DLDD1P, DLDD2P, DYDD1P, DYDD2P, DPMDD1P, DPMDD2P, 

DYMDD1P, DYMDD2P 

INPUT #1 , DRMDD1P, DRMDD2P 

INPUT #1, DLDY1P, DLDY2P, DDDY1P, DDDY2P, DPMDY1P, DPMDY2P, 

DYMDY IP , DYMDY 2 P 

INPUT #1, DRMDY1P, DRMDY2P 

INPUT # 1 , DLDPM1P , DLDPM2P, DDDPM1P , DDDPM2P , DYDPM1P , DYDP 
M2P, DYMDPM1P, DYMDPM2P 
INPUT # 1 , DRMDPM1P , DRMDPM2 P 

INPUT # 1 , DLDYM1P , DLDYM2 P , DDDYM1P , DDDYM2 P , DYDYM1P , DYDY 
M2 P , DPMDYM1P , DPMDYM2 P 

Figure A. 11 RED. BAS - Data Reduction Program 
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INPUT # 1 , DRMDYM1P , DRMDYM2 P 

INPUT # 1 , DLDRM1P , DLDRM2P, DDDRM1P , DDDRM2 P , DYDRM1P , DYDR 
M2P, DPMDRM1P, DPMDRM2P 
INPUT # 1 , DYMDRM1P , DYMDRM2 P 

INPUT #1, DDDL1N, DDDL2N, DYDL1N, DYDL2N , DPMDL1N , DPMDL2N , 

DYMDL1N, DYMDL2N 

INPUT # 1 , DRMDL1N , DRMDL2N 

INPUT # 1 , DLDD1N , DLDD2N , DYDD1N , DYDD2N, DPMDD1N , DPMDD2N , 

DYMDD1N, DYMDD2N 

INPUT # 1 , DRMDD1N , DRMDD2N 

INPUT # 1 , DLDY1N , DLDY2N, DDDY1N , DDDY2N, DPMDY1N , DPMDY2N , 

DYMDY1N, DYMDY2N 

INPUT # 1 , DRMDY IN , DRMDY 2N 

INPUT # 1 , DLDPM1N , DLDPM2N , DDDPM1N , DDDPM2N , DYDPM1N , DYDP 
M2N, DYMDPM1N 

INPUT # 1 , DYMDPM2N , DRMDPM1N , DRMDPM2N 

INPUT # 1 , DLDYM1N , DLDYM2N , DDDYM1N , DDDYM2N , DYDYM1N , DYDY 
M2N , DPMDYM1N 

INPUT # 1 , DPMDYM2N , DRMDYM1N , DRMDYM2N 

INPUT # 1 , DLDRM1N , DLDRM2N , DDDRM1N , DDDRM2N, DYDRM1N , DYDR 
M2N , DPMDRM1N 

INPUT # 1 , DPMDRM2N , DYMDRM1N , DYMDRM2N 
CLOSE #1 

» CONVERT RAW COUNTS TO FORCES 

FOR J = 1 TO NOD% 

A = ( (INCALL/ (CLL-ZL) ) * (L(J) -ZL) ) 

B = ( (INCALD/ (CALD-ZD) ) * (D(J) -ZD) ) 

C = ( (INCALY/ (CALY-ZY) ) * (Y ( J) -ZY) ) 

D=( (INCALPM/ (CALPM-ZPM) ) * (FM(J) -ZPM) ) 

E=( (INCALYM/ (CALYM-ZYM) ) * (YM(J) -ZYM) ) 

F=( ( INCALRM/ ( CALRM-ZRM) ) * (RM(J) -ZRM) ) 

IF A <= 0 THEN K1=K1LNEG:K2=K2LNEG ELSE K1=K1LP0S 
:K2=K2LPOS 

L ( J) = (K1*A) + (K2*A A 2) 

IF B <= 0 THEN K1=K1DNEG:K2=K2DNEG ELSE K1=K1DP0S 
:K2=K2DPOS 

D ( J) = (K1*B) + (K2*B A 2) 

IF C <= 0 THEN K1=K1YNEG : K2=K2 YNEG ELSE K1=K1YP0S 
: K2=K2 YPOS 

Y ( J) = (K1*C) + (K2*C A 2) 

IF D <= 0 THEN K1=K1PMNEG: K2=K2PMNEG ELSE K1=K1 
PMPOS : K2=K2PMPOS 

PM ( J) = (K1*D) + (K2*D A 2) 

IF E <= 0 THEN K1=K1 YMNEG : K2 =K2 YMNEG ELSE K1=K1 

YMPOS : K2=K2YMPOS 

YM ( J) = (K1*E) + (K2*E A 2 ) 

IF F <= 0 THEN K1=K1RMNEG: K2=K2RMNEG ELSE K1=K1 

RMPOS : K2=K2RMPOS 

RM ( J) = (K1*F) + (K2*F A 2) 

NEXT J 
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' PERFORM INTERACTON CORRECTIONS 

FOR K = 1 TO NOD% 

' FIRST CONSTANT DETERMINATION 

IF L(K) > 0 THEN GOTO 710 ELSE GOTO 720 
DD1=DDDL1P : DD6=DDDL2 P : DY1=DYDL1P : DY6=DYDL2 P : DPM1=DP 
MDL1P : DPM6=DPMDL2 P : DYM1=DYMDL1P : DYM6=DYMDL2 P : DRM1=D 
RMDL1P : DRM6=DRMDL2 P 
GOTO 730 

DD1=DDDL1N : DD6=DDDL2N : DY1=DYDL1N : DY6=DYDL2N : DPM1=DP 
MDL1N : DPM6=DPMDL2N : DYM1=DYMDL1N : DYM6=DYMDL2N : DRM1=D 
RMDL1N t DRM6=DRMDL2N 

IF D (K) >0 THEN GOTO 740 ELSE GOTO 750 
DL1=DLDD1P: DL6=DLDD2P: DY2=DYDD1P : DY7=DYDD2P: DPM2=DP 
MDD1P : DPM7=DPMDD2P: DYM2=DYMDD1P : DYM7=DYMDD2P : DRM2=D 
RMDD1P : DRM7=DRMDD2 P 
GOTO 760 

DL1=DLDD1N : DL6=DLDD2N : DY2=DYDD1N : DY7=DYDD2N : DPM2=DP 
MDD1N : DPM7=DPMDD2N : DYM2=DYMDD1N : DYM7=DYMDD2N : DRM2=D 
RMDD1N: DRM7=DRMDD2N 

IF Y (K) >0 THEN GOTO 770 ELSE GOTO 780 
DL2=DLDY1P : DL7=DLDY2P: DD2=DDDY1P : DD7=DDDY2P: DPM3=DP 
MDY1P : DPM8=DPMDY2P *. DYM3=DYMDY1P : DYM8=DYMDY2P : DRM3=D 
RMDY1P: DRM8=DRMDY2P 
GOTO 790 

DL2=DLDY1N: DL7=DLDY2N: DD2=DDDY1N: DD7=DDDY2N: DPM3=DP 
MDY1N : DPM8=DPMDY2N : DYM3=DYMDY1N : DYM8=DYMDY2N : DRM3=D 
RMDY1N: DRM8=DRMDY2N 

IF PM (K) > 0 THEN GOTO 800 ELSE GOTO 810 
DL3=DLDPM1P: DL8=DLDPM1P : DD3=DDDPM1P: DD8=DDDPM2P : DY3 
=DYDPM1P: DY8=DYDPM2P : DYM4=DYMDPM1P : DYM9=DYMDPM2P : DR 
M4=DRMDPM1P: DRM9=DRMDPM2P 
GOTO 820 

DL3=DLDPM1N:DL8=DLDPM1N:DD3=DDDPM1N:DD8=DDDPM2N:DY3 
=DYDPM1N : DY8=DYDPM2N : DYM4=DYMDPM1N : DYM9=DYMDPM2N : DR 
M4=DRMDPM1N: DRM9=DRMDPM2N 

IF YM(K) >0 THEN GOTO 830 ELSE GOTO 840 
DL4=DLDYM1P : DL9=DLDYM2P : DD4=DDDYM1P : DD9=DDDYM2P : DY4 
=DYDYM1P : DY9=DYDYM2P : DPM4=DPMDYM1P : DPM9=DPMDYM2P : DR 
M5=DRMDYM1P : DRM1 0=DRMDYM2 P 
GOTO 850 

DL4=DLDYM1N : DL9=DLDYM2N : DD4=DDDYM1N : DD9=DDDYM2N : DY4 
=DYDYM1N : DY9=DYDYM2N : DPM4=DPMDYM1N : DPM9=DPMDYM2N : DR 
M5=DRMDYM1N : DRM10=DRMDYM2N 
IF RM (K) > 0 THEN GOTO 860 ELSE GOTO 870 
DL5=DLDRM1P: DL10=DLDRM2P: DD5=DDDRM1P: DD10=DDDRM2P : D 
Y5=DYDRM1P : DY10=DYDRM2P : DPM5=DPMDRM1P : DPM10=DPMDRM2 
P : DYM5=DYMDRM1P : DYM1 0=DYMDRM2 P 
GOTO 880 

DL5=DLDRM1N : DL10=DLDRM2N : DD5=DDDRM1N : DD10=DDDRM2N : D 
Y5=DYDRM1N : DY10=DYDRM2N : DPM5=DPMDRM1N : DPM10=DPMDRM2 
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N : DYM5=DYMDKM1N : DYM10=DYMDRM2N 

880 ' FIRST INTERACTION CALCULATION 

890 LI (K) =L (K) - (DL1*D (K) ) - (DL6* (D (K) A 2 ) ) 

900 D1(K)=D(K)-(DD1*L(K) ) - (DD6* (D (K) A 2 ) ) 

950 ' SECOND CONSTANT DETERMINATION 

960 IF LI (K) > 0 THEN GOTO 970 ELSE GOTO 980 
970 DD1=DDDL1P: DD6=DDDL2 P : DY1=DYDL1P : DY6=DYDL2 P : DPM1=DP 

MDL1P : DPM6=DPMDL2P : DYM1=DYMDL1P : DYM6=DYMDL2P : DRM1=D 
RMDL1P : DRM6=DRMDL2 P 
975 GOTO 990 

980 DD1=DDDL1N: DD6=DDDL2N : DY1=DYDL1N : DY6=DYDL2N : DPM1=DP 

MDL1N : DPM6=DPMDL2N : DYM1=DYMDL1N : DYM6=DYMDL2N : DRM1=D 
RMDL1N : DRM6=DRMDL2N 

990 IF D1(K) > 0 THEN GOTO 1000 ELSE GOTO 1010 
1000 DL1=DLDD1P : DL6=DLDD2P: DY2=DYDD1P : DY7=DYDD2P : DPM2=DP 

MDD1P : DPM7=DPMDD2 P : DYM2=DYMDD1P : DYM7=DYMDD2 P : DRM2=D 
RMDD1P : DRM7=DRMDD2 P 
1005 GOTO 1020 

1010 DL1=DLDD1N : DL6=DLDD2N : DY2=DYDD1N : DY7=DYDD2N : DPM2=DP 

MDD1N : DPM7=DPMDD2N : DYM2=DYMDD1N : DYM7=DYMDD2N : DRM2=D 
RMDD1N : DRM7 =DRMDD2 N 

1020 IF Y1(K) > 0 THEN GOTO 1030 ELSE GOTO 1040 
103 0 DL2=DLDY1P : DL7=DLDY2 P : DD2=DDDY1P : DD7=DDDY2 P : DPM3 =DP 

MDY1P : DPM8=DPMDY2 P : DYM3=DYMDY1P : DYM8=DYMDY2P : DRM3=D 
RMDY1P : DRM8=DRMDY2 P 
1035 GOTO 1050 

1040 DL2=DLDY1N:DL7=DLDY2N:DD2=DDDY1N:DD7=DDDY2N:DPM3=DP 

MDY1N : DPM8=DPMDY2N : DYM3=DYMDY1N : DYM8=DYMDY2N : DRM3=D . 
RMDY1N : DRM8=DRMDY2N 

1050 IF PM1(K) > 0 THEN GOTO 1060 ELSE GOTO 1080 
1060 DL3=DLDPM1P : DL8=DLDPM1P : DD3=DDDPM1P : DD8=DDDPM2 P: DY3 

=DYDPM1P : DY8=DYDPM2P: DYM4=DYMDPM1P : DYM9=DYMDPM2 P : DR 
M4=DRMDPM1P: DRM9=DRMDPM2P 
1070 GOTO 1090 

1080 DL3=DLDPM1N : DL8=DLDPM1N : DD3=DDDPM1N : DD8=DDDPM2N : DY3 

=DYDPM1N : DY8=DYDPM2N : DYM4=DYMDPM1N : DYM9=DYMDPM2N : DR 
M4=DRMDPM1N : DRM9 =DRMDPM2 N 
1090 IF YM1(K) > 0 THEN GOTO 1100 ELSE GOTO 1110 

1100 DL4=DLDYM1P : DL9=DLDYM2P: DD4=DDDYM1P: DD9=DDDYM2P: DY4 

=DYDYM1P : DY9=DYDYM2P : DPM4=DPMDYM1P : DPM9=DPMDYM2 P : DR 
M5=DRMDYM1P: DRM10=DRMDYM2P 
1105 GOTO 1120 

1110 DL4=DLDYM1N : DL9=DLDYM2N : DD4=DDDYM1N : DD9=DDDYM2N : DY4 

=DYDYM1N : DY9=DYDYM2N : DPM4=DPMDYM1N : DPM9=DPMDYM2N : DR 
M5=DRMDYM1N: DRM10=DRMDYM2N 
1120 IF RM1(K) > 0 THEN GOTO 1130 ELSE GOTO 1140 

1130 DL5=DLDRM1P: DL10=DLDRM2P: DD5=DDDRM1P: DD10=DDDRM2P: D 

Y5=DYDRM1P: DY10=DYDRM2P: DPM5=DPMDRM1P: DPM10=DPMDRM2 
P : DYM5=DYMDRM1P : DYM1 0=DYMDRM2 P 
1135 GOTO 1150 
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1140 DL5=DLDRM1N: DL10=DLDRM2N: DD5=DDDRM1N: DD10=DDDRM2N: D 
Y5=DYDRM1N : DY10=DYDRM2N : DPM5=DPMDRM1N : DPM10=DPMDRM2 
N : DYM5=DYMDRM1N : DYM10=DYMDRM2N 
1150 ' SECOND INTERACTION CALCULATION 



1160 L2(K)=L(K)-(DL1*D1(K) ) - (DL6* (D1 (K) A 2) ) 

1170 D2 (K) =D(K) -(DD1*L1(K) ) - (DD6* (D1 (K) A 2 ) ) 

1220 '——COMPARE INTERACTION CALCULATIONS- 

1230 DFL=ABS(L1(K)-L2(K) ) : DFD=ABS (D1 (K) -D2 (K) ) 

1250 IF DFL < .0005 AND DFD < .0005 THEN GOTO 1330 

1260 IF DFL > .0005 THEN L1(K)=L2(K) 

1270 IF DFD > .0005 THEN D1(K)=D2(K) 

1320 GOTO 960 

1330 L(K)=L2(K) :D(K)=D2(K) 

1340 NEXT K 

2330 ' PERFORM WEIGHT TARE CORRECTIONS 

2335 ' — DRAG AND LIFT NEED TO BE CORRECTED FOR AXIS ORIENT 
ATION— 

2340 FOR K = 1 TO NOD% 

2350 ALPHA = AOA (K) * (3 . 141593/180) 'CONVERT TO RADIANS 
2360 DF(K) =(( (-1*D(K) ) - (W*SIN (ALPHA) )) *COS (ALPHA) ) + ((( 
-1*L (K) ) - (W*COS (ALPHA) ) ) *SIN (ALPHA) ) 

2370 LF(K)=( ( (-1*L(K) )-(W*COS (ALPHA) ) )*COS (ALPHA) ) - ((( 

-1*L(K) ) -(W*SIN( ALPHA) ) ) *SIN (ALPHA) ) 

2380 NEXT K 

2390 ' SAVE REDUCED DATA TO FILE 

2395 CLS: PRINT "DATA REDUCTION COMPLETE, SAVING DATA TO "? 
FILE$ 

2400 OPEN FILE$ FOR OUTPUT AS #2 
2410 WRITE #2 , NOD% , CONFIG$ , CONF 

2420 WRITE #2,Q,W 

2421 WRITE #2 , ZL, ZD, ZY, ZPM, ZYM, ZRM, ZAOA 

2422 WRITE #2 , CLL, CALD, CALY , CALPM, CALYM, CALRM, CALAOA 

2430 FOR J = 1 TO NOD% 

2450 WRITE #2 , LF ( J) , DF ( J) , Y ( J) , PM ( J) , YM ( J) , RM ( J) , AOA ( J) 
2460 NEXT J 

2470 ' — RETURN TO MAIN PROGRAM 

2480 COMMON PD$,FD$,DX% 

2490 CHAIN PD$+"MAIN. BAS" , 2110 

2500 ' END OF REDUCTION 
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DATA EDITOR' 



GOSUB 180 

COLOR 15,1,7: KEY OFF: CLS 
LOCATE 5,34: PRINT »*OPTIONS MENU*" 

PRINT 
PRINT 

PRINT TAB (30) "1. 

PRINT 

PRINT TAB (30) "2. 

PRINT 

PRINT TAB (30) "3. 

PRINT TAB (30) " 

PRINT 

PRINT TAB (30) "4. EXIT DATA EDITOR" 

LOCATE 20,5: INPUT "ENTER YOUR CHOICE ( 1 , 2 , 3 , 4 ) » ;REP% 
ON REP% GOSUB 390,280,2530,480 
GOTO 30 



CREATE A NEW FILE" 

EDIT EXISTING FILE" 

INCREASE NUMBER OF TEST POINTS" 
IN AN EXISTING FILE" 



' INITIALIZING CONSTANTS 1 

OPTION BASE 1: KEY OFF 

DIM LF(100) , DF ( 100 ) ,YF(100) ,PM(100) ,RM(100) 

DIM YM(100) ,AOA( 100) , COMMAND (10) 

BLANK2$ = " . " 

BLANK1S = " " 

BLANK$ = " » 

COMMAND$ ( 1) = "c": COMMANDS (2) = "r" : COMMANDS (3) 

= "d" 

COMMANDS (4) = "e": COMMANDS (5) = "s": COMMANDS (6) = 
"U": COMMANDS (7) = "q" 

RETURN 



t EDIT A FILE 

CLS: COLOR 15,1,7 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME"; FILES 
FILES = FD$+FILE$ 

GOSUB 2680 ' READ IN FILE 

PAGE% = 1: COL% = 1: SAVED% = 0: OLDXPOS% = 0: 
OLDYPOS% = 0 

GOSUB 1050 'PRINT BACKGROUND LINES 

ROW% = 3: GOSUB 1740 : GOSUB 1900 'HIGHLIGHT 

FIRST ROW AND COLUMN 

LOCATE 23,30: PRINT "LISTING": BEEP 
GOSUB 1430 'PRINT COLUMNS 1-7 AND Q VALUE 
LOCATE 23,20: PRINT BLANK2$ 

GOTO 470 



t CREATE A FILE 

COLOR 15,1,7: CLS 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME" ; FILES 
FILES = FD$ + FILES 

LOCATE 7,10: INPUT "ENTER TEST CONFIGURATION" ; CONFIG$ 
LOCATE 9,10: INPUT "ENTER THE NUMBER OF TEST POINTS 
PER RUN" ; NOD% 
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0: OLDXPOS% 



0 



PAGE% = P: COL% = Is SAVED% = 

OLDYPOS% = 0 
GOSUB 1050 'PRINT BACKGROUND LINES 
ROW% = 3: GOSUB 1740: GOSUB 1900 'HIGHLIGHT 

FIRST ROW AND COLUMN 
GOSUB 490 ' INSERT NEW DATA 

COMMON FD$ , PD$ , DX% 

CHAIN PD$+ "MAIN. BAS", 2110 
RETURN 

' CHANGE OR INSERT DATA 

LOCATE 23,12: PRINT BLANK$ 

LOCATE 23,12: INPUT REP$: IF REP$ = "" THEN GOTO 530 
GOTO 570 

LOCATE 23,12: PRINT BLANK$ : LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 
FOR I = 1 TO 500 STEP 1 
NEXT I 
GOTO 500 

REP2 $ = RIGHT$ (REP$ , 1) 

Tl% = ASC (REP2$) 

IF Tl% >= 43 AND Tl% <= 57 THEN GOTO 650 
Tl% = Tl% OR 32: REP2 $ = CHR$(T1%) 

GOSUB 2000 

IF VALID% = 0 THEN GOTO 500 ELSE VALID% = 0 
ON FLAG% GOSUB 1780,1670,2080,2340,2770,2170,2260 
IF REP2$ <> "e" THEN GOTO 500 ELSE GOTO 1040 
WHILE REP2$ <> "e" 

NEWDATA = VAL (REP$ ) 

IF QUIT% = 1 OR OLDXPOS% = 0 THEN QUIT% = 0: GOTO 
700 ELSE 

LOCATE OLDXPOS% , OLDYPOS% , 0 
PRINT USING "####.##"; OLDDATA 
LOCATE 23,12: PRINT BLANK$ 

R% = ROW% - 2 

IF COL% = 1 THEN LF (R%+18* (PAGE%-1) ) = NEWDATA: 

YPOS% = 7: GOTO 790: ELSE 

IF COL% = 2 THEN DF (R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 18: GOTO 790: ELSE 

IF COL% = 3 THEN YF (R%+18* (PAGE%-1) )=NEWDATA: YPOS% 
= 29: GOTO 790: ELSE 

IF COL% = 4 THEN PM(R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 40: GOTO 790: ELSE 

IF COL% = 5 THEN RM(R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 51: GOTO 790: ELS 

IF COL% = 6 THEN YM (R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 62: GOTO 790: ELS 

IF COL% =7 THEN AOA (R%+18* (PAGE%-1) )=NEWDATA: YPOS% 
= 71: GOTO 790: ELS 
XPOS% = ROW% 

LOCATE XPOS%,YPOS%: COLOR 0,10 

Figure A. 12 ADATA.BAS (cont.) 
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PRINT USING »####.##»; NEWDATA: MODIFIED% = 1 
SAVED% = 0: COLOR 15,1,7 

OLDYPOS% = YPOS% : OLDXPOS% = XPOS%: OLDDATA = 

NEWDATA 

FWD% = 1 

1% = R0W%+1: GOSUB 1690 'HIGHLIGHT NEW ROW 
IF 1% > NOD%+2 THEN FWD% = 1 ELSE FWD% = 0 
IF 1% > 20 THEN FWD% = 1 ELSE FWD% - 0 
GOSUB 1800 'HIGHLIGHT NEW ROW 
COLOR 15,1,7 

LOCATE 23,12: PRINT BLANK$ 

LOCATE 23, 12: INPUT REP$ : IF REP$ = "" THEN GOTO 930 
GOTO 950 

LOCATE 23,12: PRINT BLANK$ : LOCATE 23,12: PRINT 
''INVALID COMMAND": BEEP 
GOTO 900 

REP2 $ = RIGHT? (REP$,1) 

Tl% = ASC (REP2 $ ) 

IF Tl% >= 43 AND Tl% <= 57 THEN GOTO 650 
Tl% = Tl% OR 32: REP2$ = CHR$(T1%) 

GOSUB 2000 'CHECK COMMANDS 
IF VALID% = 0 THEN GOTO 500 ELSE VALID% = 0 
ON FLAG% GOSUB 1780,1670,2080,2340,2770,2170,2260 
'COMMAND SUBROUTINES 

IF REP2$ <> "e" THEN GOTO 900 ELSE GOTO 1030 
WEND 
RETURN 

PRINT BACKGROUND LINES 



COLOR 15,1, 


7 


CLS 






LOCATE 


2,1 




FOR 1% 


= 1 


TO 80 : 


LOCATE 


1,9: 


PRINT 


LOCATE 


1,20 


: PRINT 


LOCATE 


1,32 


: PRINT 


LOCATE 


1,41 


: PRINT 


LOCATE 


1,53 


: PRINT 


LOCATE 


1, 62 


: PRINT 


LOCATE 


1,74 


: PRINT 


LOCATE 


21,1 




FOR 1% 


= 1 


TO 80 : 


FOR 1% 


= 1 


TO 22 


LOCATE 


1%, 1: 



LOCATE 1% , 80 : 

NEXT 

FOR 1% = 1 TO 21 

LOCATE 1% , 14 : 
LOCATE 1% , 25 : 
LOCATE 1% , 36 : 
LOCATE 1% , 47 : 



PRINT CHR$ (220) ; : 
"LIFT" 

"DRAG" 

"YAW" 

"PITCH" 

"ROLL" 

"YAW M." 

"AOA" 

PRINT CHR$ (220) ? : 

PRINT CHR$ (222 ) 
PRINT CHR$ (222 ) 



PRINT CHR$ ( 179 ) 
PRINT CHR$ ( 179 ) 
PRINT CHR$ ( 179 ) 
PRINT CHR$ ( 179 ) 



NEXT 



NEXT 
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LOCATE 1% , 58 : PRINT CHR$(179) 

LOCATE 1% , 69 : PRINT CHR$(179) 

NEXT 

FOR 1% = 1 TO 18 
LOCATE I%+2 , 2 

PRINT USING "###" ;I%+18* (PAGE%-1) ; : PRINT 
NEXT 

COLOR 15,4 

LOCATE 22,15: PRINT "C COLUMN R ROW D DOWN E EXIT" 
LOCATE 22,46: PRINT " S SAVE U UP Q TUNNEL 
SPEED" 

COLOR 14,0: LOCATE 23,44: PRINT CONFIG$ 

COLOR 15,1,7 

LOCATE 22,2: PRINT " [Q =" : LOCATE 22,12: PRINT "]" 
LOCATE 23,3: PRINT "COMMAND :" 

RETURN 



i PRINT COLUMNS 1-7 AND Q VALUE 

FOR 1% = 1 TO 18 

LOCATE I%+2,7: PRINT USING "####.##"; LF ( I%+18 * (PAG 
E%~ 1) ) 

NEXT 

FOR 1% = 1 TO 18 

LOCATE I%+2 , 18 : PRINT USING "####.##" ;DF (I%+18* (PA 
GE%-1) ) 

NEXT 

FOR -1% = 1 TO 18 

LOCATE I%+2,29: PRINT USING "####.##" ;YF(I%+18* (PA 
GE%-1) ) 

NEXT 

FOR 1% = 1 TO 18 

LOCATE I%+2,40: PRINT USING "####.##"; PM ( I%+ 18 * (PA 
GE%-1) ) 



NEXT 

FOR 1% = 1 TO 18 

LOCATE I%+2 , 51 : PRINT USING "####.##" ;RM(I%+18* (PA 
GE%-1) ) 

NEXT 

FOR 1% = 1 TO 18 

LOCATE I%+2,62: PRINT USING "####.##" ;YM(I%+18* (PA 
GE%-1) ) 

NEXT 

FOR 1% = 1 TO 18 

LOCATE I%+2 , 71 : PRINT USING "####.##" ;AOA ( I%+18* (P 
AGE%-1) ) 

NEXT 

LOCATE 2 2, 7: PRINT USING "###";Q 
RETURN 



' ROW INDEXING 

FWD% = VAL(REP$) : IF FWD% = 0 THEN FWD% = 1 
OLDROW% = ROW% : ROW% = ROW% + FWD% 
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1700 

1710 

1720 

1730 



IF ROW% > NOD%+2 THEN ROW% = 3 
IF ROW% >20 THEN ROW% = 3 
IF ROW% < 3 THEN ROW% = 3 

LOCATE OLDROW% , 2 : PRINT USING "###" ; OLDROW%-2+18* (PAG 



i COLUMN INDEXING 

FWD% = VAL (REP$ ) : IF FWD% = 0 THEN FWD% = 1 



E%-1) ; : PRINT ":" 

1740 COLOR 14,4: LOCATE ROW%,2 

1750 PRINT USING "###" ;ROW%-2+18* (PAGE%-1) ; : PRINT " : " 

1760 COLOR 15,1,7 
1770 RETURN 
1780 
1790 

1800 OLDCOL% = COL%: COL% = COL% + FWD% 

1810 IF COL% > 7 THEN COL% = 1 
1820 IF COL% < 1 THEN COL% = 7 

1830 IF OLDCOL% = 1 THEN LOCATE 1,9: PRINT "LIFT”: GOTO 
1900 ELSE 

IF OLDCOL% = 2 THEN LOCATE 1,20: PRINT "DRAG": GOTO 
1900 ELSE 

IF OLDCOL% = 3 THEN LOCATE 1,32: PRINT "YAW": GOTO 
1900 ELSE 

IF OLDCOL% = 4 THEN LOCATE 1,41: PRINT "PITCH": GOTO 
1900 ELSE 

IF OLDCOL% = 5 THEN LOCATE 1,53: PRINT "ROLL": GOTO 
1900 ELSE 

IF OLDCOL% = 6 THEN LOCATE 1,62: PRINT "YAW M.": GOTO 
1900 ELSE 

IF OLDCOL% = 7 THEN LOCATE 1,74: PRINT "AOA" : GOTO 
1900 ELSE 
1900 COLOR 14,4 



1840 

1850 

1860 

1870 

1880 

1890 



1910 


IF COL% 
ELSE 


— 


1 


THEN 


LOCATE 


1,9: PRINT " 


LIFT": 


1920 


IF COL% 
ELSE 


= 


2 


THEN 


LOCATE 


1,20: PRINT 


"DRAG" 


1930 


IF COL% 
ELSE 


= 


3 


THEN 


LOCATE 


1,32:- PRINT 


"YAW" :' 


1940 


IF COL% 
ELSE 


= 


4 


THEN 


LOCATE 


1,41: PRINT 


"PITCH 


1950 


IF COL% 
ELSE 


= 


5 


THEN 


LOCATE 


1,53: PRINT 


"ROLL" 


1960 


IF COL% = 
1980 ELSE 


6 


THEN 


LOCATE 


1,62: PRINT 


"YAW M 


1970 


IF COL% 


= 


7 


THEN 


LOCATE 


1,74: PRINT 


"AOA" : 



1980 

1990 

2000 

2010 

2020 

2030 

2040 



ELSE 

COLOR 15,1,7 
RETURN 

• CHECK COMMANDS 

1 % = 0 

WHILE 1% < 7 AND VALID% = 0 
1% = 1% + 1 

IF REP2 $ = COMMAND$ ( 1% ) THEN FLAG% = 1%: VALID% = 1 
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2050 WEND 

2060 IF VALID% = 0 THEN GOSUB 2420 'ERROR MESSAGE 
2070 RETURN 

2080 ' SCROLL DOWN ■ 

2090 FWD% = VAL(REP$) : IF FWD% = 0 THEN FWD% = 1 
2100 PAGE% = PAGE% + FWD% 

2110 IF PAGE% < 1 THEN PAGE% = 1 

2120 GOSUB 1050 'PRINT BACKGROUND LINES 

2130 ROW% = 3: COL% = Is GOSUB 1740: GOSUB 1900 'HIGHLIGHT 
FIRST ROW AND COLUMN 

2140 GOSUB 1430 'PRINT COLUMNS 1-7 AND Q VALUE 
2150 OLDXPOS% = 0 
2160 RETURN 

217 0 ' SCROLL UP 

2180 FWD% = VAL (REP$ ) : IF FWD% = 0 THEN FWD% = 1 
2190 PAGE% = PAGE% - FWD% 

2200 IF PAGE% < 1 THEN PAGE% = 1 

2210 GOSUB 1050 'PRINT BACKGROUND LINES 

2220 ROW% = 3: COL% = 1: GOSUB 1740: GOSUB 1900 'HIGHLIGHT 
FIRST ROW AND COLUMN 

2230 GOSUB 1430 ' PRINT COLUMNS 1-7 AND Q VALUE 
2240 OLDXPOS% = 0 
2250 RETURN 

22 60 ' PRINT OUT Q VALUE 

2270 LOCATE 23,12: PRINT BLANK $ 

2280 LOCATE 23,12: INPUT "Q =";Q 

2290 LOCATE 23,12: PRINT BLANK$ 

2300 COLOR 0,10 

2310 LOCATE 22,7: PRINT USING "###"; Q 
2320 COLOR 15,1,7 
2330 RETURN 

2340 '- EXIT EDITOR 

2350 IF SAVED% = 1 OR MODIFIED% = 0 THEN GOTO 2410 

23 6 0 LOCATE 23,20: PRINT BLANK2 $ : BEEP ' 

2370 LOCATE 23,20: INPUT ''SAVE FILE (Y/N) " ;REP$ 

2380 IF REP$ = "N" OR REP$ = "n" THEN GOTO 2410 
2390 GOSUB 2770 'SAVE FILE 
2400 LOCATE 23,20: PRINT BLANK2$ 

2410 RETURN 

2420 ' ERROR MESSAGE 

243 0 LOCATE 23,20: PRINT BLANK2 $ 

2440 LOCATE 23,20: PRINT " INVALID COMMAND" 

2450 FOR F = 300 TO 500 STEP 100 
2460 SOUND F,2 

2470 SOUND 32767,2 

2480 NEXT 

2490 FOR I = 1 TO 500 STEP 1 
2500 NEXT I 

2510 LOCATE 23,20: PRINT BLANK2$ 

2520 RETURN 
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2530 ' ADD TEST POINTS TO FILE 

2540 CLS 

2550 LOCATE 3,5: PRINT"THIS OPTION ALLOWS YOU TO INCREASE 
THE NUMBER OF TEST" 

2560 LOCATE 4,5: PRINT"POINTS OF AN EXISTING FILE" 

2570 LOCATE 7,5: INPUT"WHAT IS THE NAME OF YOUR FILE" ; FILE$ 
2580 FILE$ = FD$ + FILE$ 

2590 GOSUB 2680 

2610 LOCATE 10,5: INPUT"WHAT IS THE NEW NUMBER OF TUNNEL 
RUNS" ;NOQ% 

2620 LOCATE 13,5: INPUT"WHAT IS THE NEW NUMBER OF TEST 
POINTS PER RUN" ;NOD% 

2630 GOSUB 2770 

2660 LOCATE 20,5: PRINT"NOW YOU CAN CALL UP YOUR FILE AND 
ADD IN THE NEW POINTS" 

2670 RETURN 

2 68 0 ' READ IN FILE 

2690 OPEN FILE$ FOR INPUT AS #1 
2700 INPUT # 1 , NOD% , CONFIG$ , CONF 
2710 INPUT #1,Q,W 

2714 INPUT #1 , ZL, ZD, ZY , ZPM, ZYM, ZRM, ZAOA 

2715 INPUT # 1 , CLL , CALD , CALY , CALPM , CALYM , CALRM , CALAOA 
2720 FOR J = 1 TO NOD% 

2730- INPUT # 1 , LF ( J ) , DF ( J ) , YF ( J ) , PM ( J ) , RM ( J ) , YM ( J ) , AOA ( J ) 
2740 NEXT J 
2750 CLOSE #1 
2760 RETURN 

2 770 ' SAVE FILE 

2780 LOCATE 23,20: PRINT BLANK2$ : LOCATE 23,30: PRINT 
"SAVING FILE" 

2790 OPEN FILE$ FOR OUTPUT AS #2 
2800 WRITE #2 , NOD% , CONFIG$ , CONF 
2810 WRITE #2 , Q , W 

2814 WRITE #2, ZL, ZD, ZY, ZPM, ZYM, ZRM, ZAOA 

2815 WRITE #2, CLL, CALD, CALY, CALPM, CALYM, CALRM, CALAOA 
2820 FOR J = 1 TO NOD% 

2830 WRITE #2 , LF ( J) , DF ( J) , YF ( J) , PM ( J) , RM ( J) , YM ( J) , AOA ( J) 
2840 NEXT J 
2850 CLOSE #2 
2860 RETURN 
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10 REM COMP. BAS; CALCULATE Cd, Cl AND E. F. A.— — 

15 REM 

16 DIM LF(100) , DF(IOO) , YF(IOO) , PM (100) ,YM(100) , RM(IOO) , 
AO A (100) 

17 DIM CD ( 100 ) , CL( 100) , EFA ( 100) 

20 COLOR 15,1,4: KEY OFF: CLS 

30 PRINT "WHAT IS THE NAME OF YOUR FILE CONTAINING THE"; 
35 PRINT " CALIBRATED" 

40 INPUT "WIND TUNNEL DATA" ; FILE$ 

50 FILE$ = FD$+ FILE$ 

60 GOSUB 680 'READ IN DATA FROM FILE 

70 REM — FOR THE FOLLOWING, AREAL IS MEASURED IN SQUARE 
FEET 

90 AREAL =60/144 

270 FOR J = 1 TO NOD% 

280 CL ( J) = LF ( J) / (Q*AREAL) 

290 CD ( J) = DF ( J) / (Q*AREAL) 

295 EFA ( J) =LF(J)/Q 

300 NEXT J 

31 0 ' STORE DATA FOR Cd VS. Cl 

311 INPUT "WHAT IS THE NAME FOR THE FILE TO STORE Cd VS. 

Cl DATA" ;N$ 

312 Y$="Cd":X$="Cl" 

320 OPEN FD$+N$ FOR OUTPUT AS #1. 

330 WRITE #1 , NOD% , Y$ , X$ , CONFIG$ , CONF 

340 WRITE #1 , Q 

350 FOR J = 1 TO NOD% 

360 WRITE #1 , CD ( J) , CL ( J) 

370 NEXT J 

380 CLOSE #1 

390 ' STORE DATA FOR Cd VS. Cl* Cl 

391 INPUT "WHAT IS THE NAME FOR THE FILE TO STORE Cd VS. 

Cl A 2 DATA" ;N$ 

392 Y$="Cd" :X$="C1 A 2" 

400 OPEN FD$+N$ FOR OUTPUT AS #2 
410 WRITE #2 , NOD% , Y$ , X$ , CONFIG$ , CONF 
420 WRITE #2 , Q 
430 FOR J = 1 TO NOD% 

440 B = CL ( J) *CL ( J) 

450 WRITE #2 , CD ( J) ,B 

460 NEXT J 

470 CLOSE #2 

43 0 • STORE DATA FOR Cl VS. AOA 

481 INPUT "WHAT IS THE NAME FOR THE FILE TO STORE THE Cl 

VS. AOA DATA " ; N $ 

482 Y$="C1" :X$="AOA" 

490 OPEN FD$+N$ FOR OUTPUT AS #2 
500 WRITE #2 , NOD% , Y$ , X$ , CONFIG$ , CONF 
510 WRITE #2 , Q 
520 FOR J = 1 TO NOD% 

Figure A. 13 COMP. BAS - Program to Calculate Parameters 
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WRITE #2,CL(J) ,AOA(J) 

NEXT J 
CLOSE #2 

• STORE DATA FOR E.F.A. VS AOA 

INPUT "WHAT IS THE NAME FOR THE FILE TO STORE THE 
. F. A. VS. AOA DATA " ; N $ 

Y$="E.F.A. " :X$="AOA" 

OPEN FD$+N$ FOR OUTPUT AS #1 
WRITE #1 , NOD% , Y$ , X$ , CONFIG$ , CONF 
WRITE #1 , Q 

FOR J = 1 TO NOD% 

WRITE #1,EFA(J) ,AOA(J) 

NEXT J 
CLOSE #1 

' RETURN TO MAIN PROGRAM 

COMMON FD$,PD$,DX% 

CHAIN PD$+"MAIN. BAS", 2110 

i READ IN DATA FROM CALIBRATED FILE 

OPEN FILE$ FOR INPUT AS #2 
INPUT #2 , NOD% , CONFIG$ , CONF 
INPUT #2 , Q , W 

INPUT #2 , ZL , ZD , ZY , ZPM, ZYM, ZRM, ZAOA 
INPUT # 2 , CLL , CALD , CALY , CALPM , CALYM , CALRM , CALAOA 
FOR J = 1 TO NOD% 

INPUT #2 , LF ( J) , DF ( J) ,YF(J) ,PM(J) ,RM(J) ,YM(J) , 
AOA(J) 

NEXT J 
CLOSE #2 
RETURN 



Figure A. 13 COMP. BAS (cont.) 
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' DATA EDITOR 

GOSUB 180 

COLOR 15,1,7: KEY OFF: CLS 

LOCATE 5,34: PRINT "*OPTIONS MENU*" 

PRINT 

PRINT 



PRINT TAB (30) 


"1. 


CREATE A NEW FILE" 


PRINT 






PRINT TAB (30) 


"2 . 


EDIT EXISTING FILE" 


PRINT 






PRINT TAB (30) 


"3. 


INCREASE NUMBER OF TEST POINTS" 


PRINT TAB (30) 


ii 


IN AN EXISTING FILE" 


PRINT 






PRINT TAB (30) 


"4. 


EXIT DATA EDITOR" 


LOCATE 20,5: 


INPUT 


"ENTER YOUR CHOICE ( 1 , 2 , 3 , 4 ) " ;REP% 



ON REP% GOSUB 380,270,2200,490 
GOTO 30 

' INITIALIZING CONSTANTS 

OPTION BASE 1: KEY OFF 

DIM Y ( 100) , X ( 100) , COMMAND ( 10) 

BLANK2 $ = " "to 

BLANK 1$ = " " 

BLANK$ = " » 

COMMAND$ ( 1) = "c": COMMAND? (2) = "r": COMMAND? (3) = "d" 
COMMAND? (4) = "e": COMMAND? (5) = »s": COMMAND? (6) = 
"U": COMMAND? (7) = "q" 

RETURN 

' — EDIT A FILE 

CLS: COLOR 15,1,7 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME" ; FILE? 

FILE? = FD? + FILE? 

GOSUB 2350 ' READ IN FILE 

PAGE% = 1: COL% =1: SAVED% = 0: OLDXPOS% = 0: 

OLDYPOS% = 0 

GOSUB 1010 'PRINT BACKGROUND LINES 

ROW% = 3: GOSUB 1520 : GOSUB 1620 'HIGHLIGHT 

FIRST ROW AND COLUMN 

LOCATE 23,30: PRINT "LISTING": BEEP 

GOSUB 1360 'PRINT COLUMNS 1-7 AND Q VALUE 

LOCATE 23,20: PRINT BLANK2? 

GOTO 480 

' CREATE A FILE 

COLOR 15,1,7: CLS 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME"; FILE? 

FILR? = FD? + FILE? 

LOCATE 7,10: INPUT "ENTER TEST CONFIGURATION" ; CONFIG? 
LOCATE 9,10: INPUT "ENTER THE NUMBER OF TEST POINTS 
PER RUN" ;NOD% 

LOCATE 11,10: INPUT "ENTER TITLE FOR X-VALUES (i.e. 

CL, AOA) " ;X1? 

Figure A. 14 BDATA.BAS - On-Screen Data Editor 
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LOCATE 13,10: INPUT "ENTER TITLE FOR Y-VALUES (i.e. 

Cd, Cl) " ; ¥1$ 

PAGE% = 1: COL% = 1: SAVED% = 0: OLDXPOS% = 0: 

OLDYPOS% = 0 

GOSUB 1010 'PRINT BACKGROUND LINES 

ROW% = 3: GOSUB 1520: GOSUB 1620 'HIGHLIGHT FIRST 

ROW AND COLUMN 

GOSUB 500 ' INSERT NEW DATA 

COMMON FD$ , PD$ , DX% 

CHAIN PD$+"MAIN. BAS", 2110 
RETURN 

i CHANGE OR INSERT DATA 

LOCATE 23,12: PRINT BLANK$ 

LOCATE 23,12: INPUT REP$: IF REP$ = "" THEN GOTO 540 
GOTO 580 

LOCATE 23,12: PRINT BLANK$ : LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 
FOR I = 1 TO 500 STEP 1 
NEXT I 
GOTO 510 

REP 2$ = RIGHT$ (REP$ , 1) 

Tl% = ASC (REP2$) 

IF Tl% >= 43 AND Tl% <= 57 THEN GOTO 660 
Tl% = Tl% OR 32: REP2$ = CHR$(T1%) 

GOSUB 1670 

IF VALID% = 0 THEN GOTO 510 ELSE VALID% = 0 
ON FLAG% GOSUB 1550,1450,1750,2010,2440,1840,1930 
IF REP2$ <> "e" THEN GOTO 510 ELSE GOTO 1000 
WHILE REP2$ <> "e" 

NEWDATA = VAL (REP$) 

IF QUIT% = 1 OR OLDXPOS% = 0 THEN QUIT% = 0: GOTO 
710 ELSE 

LOCATE OLDXPOS% , OLDYPOS% , 0 
PRINT USING "####.##",* OLDDATA 
LOCATE 23,12: PRINT BLANK$ 

R% = ROW% - 2 

IF COL% = 1 THEN Y (R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 18: GOTO 750: ELS 

IF COL% - 2 THEN X (R%+18* (PAGE%-1) ) =NEWDATA: YPOS% 
= 43: GOTO 750: ELSE 
XPOS% = ROW% 

LOCATE XPOS% , YPOS% : COLOR 0,10 

PRINT USING "####.##"; NEWDATA: MODIFIED% = 1 
SAVED% = 0: COLOR 15,1 

OLDYPOS% = YPOS% : OLDXPOS% = XPOS% : OLDDATA = 
NEWDATA 
FWD% = 1 

1% « R0W%+1 : GOSUB 1470 'HIGHLIGHT NEW ROW 
IF 1% > NOD%+2 THEN FWD% = 1 ELSE FWD% = 0 
IF 1% > 20 THEN FWD% = 1 ELSE FWD% = 0 

Figure A. 14 BDATA.BAS (cont.) 
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840 GOSUB 1570 'HIGHLIGHT NEW ROW 
850 COLOR 15,1,7 

860 LOCATE 23,12: PRINT BLANK? 

870 LOCATE 23,12: INPUT REP$: IF REP$ = "" THEN GOTO 

890 

880 GOTO 910 

890 LOCATE 23,12: PRINT BLANK?: LOCATE 23,12: PRINT 

"INVALID COMMAND": BEEP 
900 GOTO 860 

910 REP2? = RIGHT? (REP$,1) 

920 Tl% = ASC (REP2?) 

930 IF Tl% >= 43 AND Tl% <= 57 THEN GOTO 660 

940 Tl% = Tl% OR 32: REP2? = CHR?(T1%) 

950 GOSUB i.670 'CHECK COMMANDS 

960 IF VALID% = 0 THEN GOTO 510 ELSE VALID% = 0 

970 ON FLAG% GOSUB 1550,1450,1750,2010,2440,1840,1930 

'COMMAND SUBROUTINES 

980 IF REP 2? <> "e" THEN GOTO 860 ELSE GOTO 990 

990 WEND 

1000 RETURN 

1010 ' — PRINT BACKGROUND LINES 

1020 COLOR 15,1,7 

1030 CLS 

1040 LOCATE 2,1 

1050 FOR 1% = 1 TO 80 : PRINT CHR?(220);: NEXT 
1060 LOCATE 1,18: PRINT »Y = "?Y1$ 

1070 LOCATE 1,43: PRINT "X = ";X1? 

1080 LOCATE 1,68: PRINT "COMMANDS" 

1090 LOCATE 21,1 

1100 FOR 1% = 1 TO 80 : PRINT CHR?(220);: NEXT 
1110 FOR 1% = 1 TO 22 

1120 LOCATE 1% , 1 : PRINT CHR?(222) 



1130 


LOCATE 


1%, 80: 


PRINT 


CHR? ( 222 ) 


1140 


NEXT 








1150 


FOR 1% = 1 


TO 21 






1160 


LOCATE 


1% , 28 : 


PRINT 


CHR? (179) 


1170 


LOCATE 


1% , 53 : 


PRINT 


CHR? (179) 


1180 


NEXT 








1190 


FOR 1% = 1 


TO 18 






1200 


LOCATE 


I%+2,2 







1210 PRINT USING "###"; I%+18 * (PAGE%-1) ; : PRINT ":" 

1220 NEXT 



1230 


COLOR ; 


15,4 










1240 


LOCATE 


3 , 66 


PRINT 


"C 


COLUMN 


ii 


1250 


LOCATE 


4,66 


PRINT 


"R 


ROW 




1260 


LOCATE 


5,66 


PRINT 


»D 


DOWN 




1270 


LOCATE 


6,66 


PRINT 


"E 


EXIT 




1280 


LOCATE 


7,66 


PRINT 


"S 


SAVE 




1290 


LOCATE 


8,66 


PRINT 


"U 


UP 




1300 


LOCATE 


9,66 


PRINT 


"Q 


TUNNEL 


SPEED" 



Figure A. 14 BDATA.BAS (cont.) 
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1310 COLOR 14,0: LOCATE 23,44: PRINT CONFIG$ 

1320 COLOR 15,1,7 

1330 LOCATE 22,2: PRINT "[Q =" : LOCATE 22,12: PRINT "]" 

1340 LOCATE 23,3: PRINT "COMMAND :" 

1350 RETURN 

1360 ' PRINT COLUMNS 1-7 AND Q VALUE 

1370 FOR 1% = 1 TO 18 

1380 LOCATE I%+2,18: PRINT USING »####.##" ;Y (I%+18* (PAG 
E%-1) ) 

1390 NEXT 

1400 FOR 1% = 1 TO 18 

1410 LOCATE I%+2 , 43 : PRINT USING "####.##"; X ( I%+18* (PAG 

E%-1) ) 

1420 NEXT 

143 0 LOCATE 2 2, 7 -.PRINT USING »###»;Q 
1440 RETURN 

1450 ' ROW INDEXING 

1460' FWD% = VAL (REP$ ) : IF FWD% = 0 THEN FWD% =’ 1 
1470 OLDROW% = ROW% : ROW% = ROW% + FWD% 

1480 IF ROW% > NOD%+2 THEN ROW% = 3 
1490 IF ROW% >20 THEN ROW% = 3 
1500 IF ROW% < 3 THEN ROW% = 3 

1510 LOCATE OLDROW% , 2 : PRINT USING "###"; OLDROW%-2+18 *( PAG 
E%-1) ; : PRINT 

1520 COLOR 14,4: LOCATE ROW% , 2 : PRINT USING "###"; ROW%-2+l 
8* (PAGE%-1) ; : PRINT " : " 

1530 COLOR 15,1,7 
1540’ RETURN 

1550 ' COLUMN INDEXING 

1560 FWD% = VAL(REP$) : IF FWD% = 0 THEN FWD% = 1 
1570 OLDCOL% = COL% : COL% = COL% + FWD% 

1580 IF COL% > 2 THEN COL% = 1 
1590 IF COL% < 1 THEN COL% = 2 

1600 IF OLDCOL% = 1 THEN LOCATE 1,18: PRINT "Y = "?Y1$: 

GOTO 1620 ELSE 

1610 IF OLDCOL% = 2 THEN LOCATE 1,43: PRINT "X = ";X1$: 

GOTO 1620 ELSE 
1620 COLOR 14,4 

1630 IF COL% = 1 THEN LOCATE 1,18: PRINT "Y = "?Y1$: GOTO 
1650 ELSE 

1640 IF COL% = 2 THEN LOCATE 1,43: PRINT "X = "?X1$: GOTO 
1650 ELSE 
1650 COLOR 15,1,7 
1660 RETURN 

1670 ' CHECK COMMANDS 

1680 1% = 0 

1690 WHILE 1% < 7 AND VALID% = 0 
1700 1% = 1% + 1 

1710 IF REP2 $ = COMMAND$ (1%) THEN FLAG% = 1%: VALID% = 1 
1720 WEND 



Figure A. 14 BDATA.BAS (cont.) 



71 



1730 IF VALID% = 0 THEN GOSUB 2090 'ERROR MESSAGE 
1740 RETURN 

1750 ' SCROLL DOWN 

1760 FWD% = VAL(REP$): IF FWD% = 0 THEN FWD% = 1 
1770 PAGE% = PAGE% + FWD% 

1780 IF PAGE% < 1 THEN PAGE% = 1 

1790 GOSUB 1010 'PRINT BACKGROUND LINES 

1800 ROW% = 3: COL% = 1: GOSUB 1520: GOSUB 1620 'HIGHLIGHT 
FIRST ROW AND COLUMN 

1810 GOSUB 1360 'PRINT COLUMNS 1-7 AND Q VALUE 
1820 OLDXPOS% = 0 
1830 RETURN 

1840 ' SCROLL UP 

1850 FWD% = VAL (REP$ ) : IF FWD% = 0 THEN FWD% = 1 
1860 PAGE% = PAGE% - FWD% 

1870 IF PAGE% < 1 THEN PAGE% = 1 

1880 GOSUB 1010 'PRINT BACKGROUND LINES 

1890 ROW% = 3: COL% = 1: GOSUB 1520: GOSUB 1620 'HIGHLIGHT 
FIRST ROW AND COLUMN 

1900 GOSUB 1360 ' PRINT COLUMNS 1-7 AND Q VALUE 
1910 OLDXPOS% = 0 
1920 RETURN 

1930 ' — PRINT OUT Q VALUE — — 

1940 LOCATE 23,12: PRINT BLANK$ 

1950 LOCATE 23,12: INPUT "Q =";Q 
1960 LOCATE 23,12: PRINT BLANK$ 

1970 COLOR 0,10 

1980 LOCATE 22,7: PRINT USING "###"’; Q 
1990 COLOR 15,1,7 
2000 RETURN 

2010 ' EXIT EDITOR 

2020 IF SAVED% = 1 OR MODIFIED% = 0 THEN GOTO 2080 

203 0 LOCATE 23,20: PRINT BLANK2 $ : BEEP 

2040 LOCATE 23,20: INPUT "SAVE FILE (Y/N)"?REP$ 

2050 IF REP$ = "N" OR REP$ = "n" THEN GOTO 2080 
2060 GOSUB 2440 'SAVE FILE 
2070 LOCATE 23,20: PRINT BLANK2$ 

2080 RETURN 

2090 '— ERROR MESSAGE 

2100 LOCATE 23,20: PRINT BLANK2$ 

2110 LOCATE 23,20: PRINT " INVALID COMMAND" 

2120 FOR F = 300 TO 500 STEP 100 
2130 SOUND F , 2 

2140 SOUND 32767,2 

2150 NEXT 

2160 FOR I = 1 TO 500 STEP 1 
2170 NEXT I 

218 0 LOCATE 23,20: PRINT BLANK2 $ 

2190 RETURN 

2200 ' ADD TEST POINTS TO FILE 

Figure A. 14 BDATA.BAS (cont.) 
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CLS 

LOCATE 3,5: PRINT"THIS OPTION ALLOWS YOU TO INCREASE 
THE NUMBER OF TEST" 

LOCATE 4,5: PRINT"POINTS OF AN EXISTING FILE" 

LOCATE 7,5: INPUT"WHAT IS THE NAME OF YOUR FILE" ? FILE$ 
FILE$ = FD$ + FILE$ 

GOSUB 2350 

LOCATE 13,5: INPUT"WHAT IS THE NEW NUMBER OF TEST 
POINTS PER RUN" ;NOD% 

GOSUB 2440 

LOCATE 20,5: PRINT"NOW YOU CAN CALL UP YOUR FILE AND 
ADD IN THE NEW POINTS" 

RETURN 

i READ IN FILE 

OPEN FILE? FOR INPUT AS #1 

INPUT #1 , NOD% , Yl$ , Xl$ , CONFIG? , CONF 

INPUT #1 , Q 

FOR J = 1 TO NOD% 

INPUT #1,Y(J) ,X(J) 

NEXT J 
CLOSE #1 
RETURN 

t SAVE FILE 

LOCATE 2 3,20: PRINT BLANK2 ? : LOCATE 23,30: PRINT 
"SAVING FILE" 

OPEN FILE? FOR OUTPUT AS #2 
WRITE #2 , NOD% , Yl? , XI? , CONFIG? , CONF 
WRITE #2 , Q 

FOR J = 1 TO NOD% 

WRITE #2,Y(J) ,X(J) 

NEXT J 
CLOSE #2 
RETURN 
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10 ' PLOTTING ROUTINE — 

20 DIM X(300) , Y ( 3 00 ) ,B(300) ,C(300) , D ( 3 00 ) 

30 KEY OFF: SCREEN 2 

40 HDPOS% =44: B$ = " " 

60 GOSUB 260 
70 OPT$ = CURVES (1) 

80 ERASE X,Y,FILESTK$, CURVES, B,C,D 

90 CHAIN PD$+"MAIN.BAS", 2110, ALL 

150 ' — READ IN FILE — — 

160 OPEN FILES FOR INPUT AS #2 
170 INPUT #2 , NOD% , Yl$ , Xl$ , H$ , CONF 
180 INPUT #2 , Q 
190 FOR 1% = 1 TO NOD% 

200 INPUT #2, Y(I%) ,X(I%) 

210 NEXT 

220 CLOSE #2 
230 RETURN 

260 ' INITIALIZING PLOTTER 

270 FOR L% = 1 TO NOF% 

280 FILES = FILESTK$ (L%) 

290 GOSUB 160 

300 IF L%=1 THEN XMIN=X(1)-: XMAX=XMIN: YMIN=Y(1): 

YMAX=YMIN 

310 FOR I%=1 TO NOD% 

320 IF X (1%) <XMIN THEN XMIN=X(I%) ELSE IF X(I%)> 

XMAX THEN XMAX=X(I%) 

330 IF Y (1%) <YMIN THEN YMIN=Y(I%) ELSE IF Y(I%)> 

YMAX THEN YMAX=Y(I%) 

340 NEXT 1% 

350 NEXT L% 

351 CLS 

352 PRINT "AFTER CURVES ARE PLOTTED YOU HAVE THE OPTION 
OF OBTAINING A HARD" 

353 PRINT "COPY BY PRESSING CTRL-PRTSC. THE CURVE WILL 
REMAIN ON THE SCREEN" 

354 PRINT "UNTIL YOU PRESS THE F2 KEY." 

355 PRINT 

360 INPUT "HEADDING FOR PLOT =";HEAD$ 

370 T1%=LEN (HEADS) : IF Tl%>60 THEN HEAD$=LEFT$ (HEADS , 60) 
380 INPUT "X TITLE FOR THE PLOT =";XTITLE$ 

390 T1%=LEN (XTITLES) : IF Tl%>60 THEN XTITLE$=LEFT$ (XTITLE 
$,60) 

400 INPUT "Y TITLE FOR THE PLOT =";YTITLE$ 

410 T1%=LEN (YTITLES) : IF Tl%>60 THEN YTITLE$=LEFT$ (YTITLE 
$,60) 

420 PRINT "MINIMUM X VALUE = [ DEFAULT XMIN ?»]"; : INPUT 
Tl$ 

430 IF Tl$ <> "" THEN TEMP=VAL ( T 1 $ ) ELSE TEMP=XMIN 

440 IF TEMP<=XMIN THEN 0K%=1 ELSE OK%=0 

450 IF Tl$ <> "" THEN TEMP3 =XMIN: XMIN = TEMP 

Figure A. 15 PLOTTER. BAS - Screen Graphics Program 
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PRINT "MAXIMUM X VALUE = [DEFAULT: " ;XMAX; "]"? : INPUT 
Tl$ 

IF Tl$ <> "" THEN XMAX= VAL ( T 1 $ ) 

PRINT "MINIMUM Y VALUE = [ DEFAULT YMIN : INPUT 

Tl$ 

IF Tl$ <> "" THEN YMIN=VAL(T1$) 

PRINT "MAXIMUM Y VALUE = [ DEFAULT YMAX : INPUT 
Tl$ 

IF Tl$ <> "" THEN YMAX=VAL ( T 1 $ ) 

RANGE=XMAX-XMIN 

IF 0K%=1 THEN PTR%=1: STKPTR%=1: PGSTK%(1)=1 ELSE 

GOSUB 580 

CLS : SCREEN 2 

GOSUB 840 

RETURN 

i SET PAGES 

TEMP 2= (TEMP -TEMP 3 ) /RANGE+1 : PG%=INT (TEMP 2 ) 

STKPTR%=PG%+ 1 

FOR 1% = 1 TO STKPTR% 

T2=TEMP-RANGE* (I%-1) 

J%=1 

WHILE X ( J%) <T2 
J%=J%+1 
WEND 

PGSTK% ( PG%-I%+2 ) =J% 

NEXT 1% 

PTR%=PGSTK% ( STKPTR% ) 

RETURN 

i p L oT WITHOUT GRID 

LINE (ORIGINX% , ORIGINY%) - (ORIGINX%+7*INTERVALX% , TEMP4% 

) ,1,B 

X%=ORIGINX% 

FOR 1% = 1 TO 8 

LINE (X%, 162) -(X%, 165) 

X%=X%+INTERVALX% 

NEXT 

Tl%=ORIGINX%-8 : T2%=ORIGINY% 

FOR 1% = 1 TO 5 

LINE (T1%,T2%) - (ORIGINX% , T2 % ) 

T2 %=T2 %+ INTERVALY % 

NEXT 

RETURN 

i PLOTTING ROUTINES 

ORIGINY%=12 : ORIGINX%=72: INTERVALY%=3 0 : INTERVALX%=80 
TEMP4%=ORIGINY%+5*INTERVALY% 

VIEW (0,0) -(639,163) : CLS: VIEW 
GOSUB 710 

IF LEFT$ (GRD$ , 1) ="G" THEN GOSUB 2400 
X%=ORIGINX% : Y%=ORIGINY% 

FOR L%=1 TO NOF% 

Figure A. 15 PLOTTER. BAS (cont.) 
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IF N0F%=1 THEN GOTO 960 
FILE$=FILESTK$ (L%) 

GOSUB 160 
XINDEX=XMIN 
GOSUB 1410 

IF L%=1 THEN GOSUB 1120 

IF L% <> NOF% THEN PTR%=PGSTK% (STKPTR% ) 

NEXT L% 

ON KEY ( 1 ) GOSUB 1006 

ON KEY (2) GOSUB 1009 

KEY ( 2 ) OFF: KEY ( 1) ON: KEY ( 2 ) ON 

IF FLAG = 1 THEN GOTO 1100 

GOTO 1001 

LOCATE 23,1: PRINT "HI" 

RETURN 
FLAG = 1 
RETURN 
FLAG = 0 
RETURN 

i p L0T ZER0 LINE 

IF (YMAX*YMIN <0) THEN TEMP= (-1*YMIN/ (YMAX-YMIN) ) *150 
ELSE GOTO 1150 
DELTAY%=CINT (TEMP) 

LINE (ORIGINX% ,TEMP4%-DELTAY%) - (ORIGINX%+7*INTERVALX% , 
TEMP4%-DELTAY%) 

IF (XMAX*XMIN <0) THEN TEMP= ( -1*XMIN/ (XMAX-XMIN) ) *550 
ELSE GOTO 1190 
DELTAX%=CINT (TEMP) 

LINE (ORIGINX%+DELTAX% , ORIGINY%) - (ORIGINX%+DELTAX% , TEM 
P4%) 

i PRINT INDECIES 

YINDEX=YMAX 

TEMP = (YMAX-YMIN) /5 

P%=LEN (HEAD$) : P%=HDPOS%- (P%* . 5) 

LOCATE 1,1: PRINT TAB(P%) HEAD$ 

FOR 1% = 1 TO 6 



PRINT USING "###.###",* YINDEX 



IF (I%= 3) 
GOTO 1260 


THEN 


FOR 


J% = 1 


TO 


2 : 


PRINT : 


NEXT 


J% 


IF ( I%<>6 ) 


THEN 


FOR 


J% = 1 


TO 


3 : 


PRINT : 


NEXT 


J% 



YINDEX=YINDEX-TEMP 
NEXT 1% 

XINDEX = XMIN: XINCR= (XMAX-XMIN) /7 
LOCATE 22,1 

PRINT USING "########.#" ;XINDEX; 

FOR 1% = 1 TO 7 

XINDEX = XINDEX + XINCR 

PRINT USING "########.#"; XINDEX; : NEXT 1% 
P%=LEN(XTITLE$) : P%=HDPOS%- (P%* . 5 ) 

LOCATE 23,1: PRINT TAB(P%) XTITLE$ 

Figure A. 15 PLOTTER. BAS (cont.) 
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T1%=LEN(YTITLE$) : P%=11- (Tl%* . 5) 

FOR 1% = 1 TO Tl% 

T1$=MID$ ( YTITLE$ , 1% , 1) : LOCATE P%+I%,1: PRINT Tl$ 
NEXT 
RETURN 

i p L oT POINTS ON GRAPH 

TEMP2=1/ (YMAX-YMIN) : TEMP3=1/ (XMAX-XMIN) 

IF MRK$=" " THEN MRK$="MARK " 

REP2 $=LEFT$ (MRK$ , 1) 

T%=L% : REP$=LEFT$(CURVE$(L%) ,1) 

IF T%>3 THEN T%=T%-3 : GOTO 1450 
IF REP$="C" THEN GOSUB 1870 
COUNT%=0: NDP%=0 

WHILE (PTR%<=NOD%) AND (X(PTR%) <= XMAX) 

IF X (PTR%) < XMIN THEN GOTO 1650 
OLDY%=Y% : OLDX%=X% 

NDP% = NDP%+1 

TEMP1=ABS ( (X (PTR%) -XMIN) *TEMP3 ) *550 
TEMP=ABS ( (Y (PTR%) -YMIN) *TEMP2) *150 
DELTAY%=CINT (TEMP) : DELTAX%=CINT (TEMPI) 

Y%=TEMP4 % -DELTAY % : X%=ORIGINX%+DELTAX% 

IF COUNT%=0 THEN OLDY%=Y%: OLDX%=X%: PSET(X%,Y%): 
FIRSTX%=X% : FIRSTY%=Y%: C0UNT%=1 
IF REP$="P" THEN 1600 
LINE- (X% , Y%) , 1 

IF REP2$ <> "M" AND REP <> "P" THEN GOTO 1630 
IF T%=1 THEN CIRCLE(X%,Y%) ,2,1,0,6.28,5/12 
IF T%=2 THEN LINE (X% , Y%) - (X%+6 , Y%) : LINE- (X%+3 , Y% 
-3): LINE- (X% , Y%) 

ELSE LINE (X% , Y% ) - (X%+6 , Y%-3 ) ,1,B : PSET(X%,Y%) 
GOTO 1640 
PSET (X% , Y%) 

IF REP$="C" THEN GOSUB 2140 
PTR%=PTR%+1 
WEND 
RETURN 

i CUBIC INTERPOLATION 

DIM H(200) , A ( 2 00) ,L(200) ,U(200) ,Z(200) 

FOR 1% = 1 TO N0D%-1 

H(I%)=X(I%+1) -X(I%) 

NEXT 

FOR I%=2 TO N0D%-1 

T1=Y(I%+1) *H(I%-1) 

T2=Y(I%) * (X(I%+1) -X(I%-1) ) 

T3=Y (I%-1) *H(I%) 

A (1%) =3 * (T1-T2+T3 ) / (H (I%-1) *H(I%) ) 

NEXT 1% 

L(l) =1: U (1) =0 : Z(l)=o 
FOR I%=2 TO NOD% 

L ( 1% ) =2 *(X(I%+1) -X(I%-1) ) -H (I%-1) *U(I%-1) 

Figure A. 15 PLOTTER. BAS (cont.) 
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U(I%)=H(I%)/L(I%) 

Z(I%)=(A(I%) -H ( I % - 1 ) *Z(I%-1) )/L(I%) 

NEXT 1% 

L(NOD%) =1: Z (NOD%) =0 : C(NOD%)=0 
FOR J%=N0D%-1 TO 1 STEP -1 

C(J%)=Z (J%) -U(J%) *C(J%+1) 

T1=(Y(J%+1)-Y(J%) )/H(J%) 

T2=H(J%) * (C(J%+1)+2*C(J%) )/3 
B(J%)=T1-T2 

D ( J% ) = (C ( J%+1) -C(J%) )/(3*H(J%) ) 

NEXT J% 

ERASE H, A, L,U, Z 
RETURN 

' PLOT CUBIC 

IF PTR%+1 > NOD% OR X(PTR%+1) > XMAX THEN GOTO 2280 
Tl= (XMAX-XMIN) /560 : T2=X(PTR%) 

TEMP1=ABS ( (X(PTR%+1) -XMIN) *TEMP3) *550 
DELTAX%=CINT (TEMPI) 

XNEXT%=ORIGINX%+DELTAX% 

X%=X%+2 

IF X% >= XNEXT% THEN GOTO 2280 

T2=T2 + 2 *T1 

T4=T2-X(PTR%) 

T3=Y ( PTR% ) +B ( PTR% ) *T4+C(PTR%) *T4*T4+D(PTR%) *T4*T4*T4 
TEMP=ABS ( (T3-YMIN) *TEMP2) *150 
DELTAY%=CINT (TEMP) : Y%=TEMP4%-DE LTA Y % 

IF Y%>17 2 THEN Y%=172 ELSE IF Y%<12 THEN Y%=12 

LINE- (X% , Y%) : GOTO 2190 

RETURN 

' PLOT WITH GRIDS 

FOR 1% = 1 TO 9 

Tl%=ORIGINY%+I%*15 

LINE (73 , Tl%) - (623 , Tl%) , , ,&H4444 

NEXT 1% 

FOR I%= 1 TO 13 

T1%=ORIGINX%+I%*40 

LINE (Tl% , 12 ) - (Tl% , 162 ) , , , &HAAAA 

NEXT 1% 

RETURN 

END 



Figure A. 15 PLOTTER. BAS (cont.) 
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NAME: Data Acquisition And Control ( DAAC) 

HEADER for BASICA 

FILE NAME: DACHDR . BAS 

DOS DEVICE NAME: DAAC 

RESERVED FUNCTION NAMES: 

AINM, AINS, AINSC, AOUM, AOUS , 

BINM, BINS, BITINS, BITOUS , BOUM, BOUS , 
CINM, CINS, CSET, DELAY 
RESERVED DEF SEG VALUE NAME: DSEG 

NAMES DEFINED AND USED BY HEADER: 

ADAPT % , AI, COUNT, FOUND%, 

HNAME$ , SG% , STAT% 



When using the BASICA Interpreter, this header 
must be executed before any function calls are 
made that access the DAAC adapter. It initializes 
a number of variables for each function call. These 
variables are reserved and should not be used except 
to access the DAAC adapter. This routine also does a 
DEF SEG to the segment where the DAAC Device Driver 
(DAC.COM) is loaded. If you execute a DEF SEG to 
access other hardware, you must DEF SEG to the segment 
of the DAAC Device Driver before any subsequent 
calls to access the DAAC adapter. 



FOUND% =0 
SG% = &H2E 

'Start searching the interrupt vectors until you find 
'one that points to the DAAC device driver. 

'Do a DEF SEG to that segment. 

WHILE ( ( SG% <= &H3E) AND (FOUND% = 0) ) 

DEF SEG = 0 

DSEG = PEEK(SG%) + PEEK(SG% + 1) * 256 

DEF SEG = DSEG 
HNAME$=" " 

FOR AI=10 TO 17 

HNAME$ = HNAME$ + CHR$ (PEEK (AI) ) 

NEXT AI 

IF HNAME$ = "DAAC " AND PEEK (18) + PEEK (19) <> 
0 THEN FOUND% = 1 
SG% = SG% + 4 



WEND 

IF FOUND% = 
NOT FOUND" : 



0 THEN PRINT "ERROR: 
END 



DEVICE DRIVER DAC.COM 



Figure A. 16 BALCAL.BAS - Balance Calibration Program 
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490 'Now initialize all function name variables for calls 
500 'to access the device driver. 



510 


AINM 


=5 


PEEK ( &H13 ) 


* 


256 


+ 


PEEK ( &H12 ) 


520 


AINS 


= 


PEEK ( &H15 ) 


* 


256 


+ 


PEEK ( &H14 ) 


530 


AINSC 


= 


PEEK(&H17) 


* 


256 


+ 


PEEK ( &H16 ) 


540 


AOUM 


= 


PEEK(&H19) 


* 


256 


+ 


PEEK ( &H18 ) 


550 


AOUS 


= 


PEEK ( &H1B) 


* 


256 


+ 


PEEK ( &H1A) 


560 


BINM 


= 


PEEK ( &H1D) 


* 


256 


+ 


PEEK ( &H1C) 


570 


BINS 


ss 


PEEK(SHIF) 


* 


256 


+ 


PEEK ( &H1E) 


580 


BITINS 


= 


PEEK ( &H2 1) 


* 


256 


+ 


PEEK(&H20) 


590 


BITOUS 


= 


PEEK ( &H23 ) 


* 


256 


+ 


PEEK ( &H22 ) 


600 


BOUM 


= 


PEEK(&H25) 


* 


256 


+ 


PEEK ( &H24 ) 


610 


BOUS 


= 


PEEK ( &H27 ) 


* 


256 


+ 


PEEK ( &H2 6 ) 


620 


CINM 


= 


PEEK ( &H29 ) 


* 


256 


+ 


PEEK ( &H28 ) 


630 


CINS 


= 


PEEK ( &H2B) 


* 


256 


+ 


PEEK ( &H2A) 


640 


CSET 


SB 


PEEK ( &H2D) 


* 


256 


+ 


PEEK ( &H2C) 


650 


DELAY 


=s 


PEEK ( &H2F) 


* 


256 


+ 


PEEK ( &H2E) 


660 


' Finally, 


execute any call 


to : 


re- 


-initialize the 



670 'device driver from any former invocation of BASIC. 

680 ADAPT % = 0 
690 COUNT = 1 
700 STAT% = 0 

710 CALL DELAY ( ADAPT %, COUNT, STAT%) 

720 ' 

730 'End of DAAC BASICA Header 
740 ' 

750 REM -BALCAL. BAS (CALIBRATE BALANCE) 

760 REM 

765 DIM L(100) ,D(100) ,Y(100) ,PM(100) ,YM(100) ,RM(100) , LBS ( 
100 ) 

769 DIM DAT (3 99) ,DAT%(399) ,DAT1(399) ,DAT1%(399) 

770 FOR I =. 1 TO 12 

78 0 IF I - 1 THEN FILE$' = "LIFTP" 

790 IF I = 2 THEN FILE? = " DRAGP" 

800 IF I = 3 THEN FILE$ = "YAWP" 

810 IF I = 4 THEN FILE? = "PITCHMP" 

820 IF I = 5 THEN FILE? = "YAWMP" 

830 IF I = 6 THEN FILE? = "ROLLMP" 

840 IF I = 7 THEN FILE? = "LIFTN" 

850 IF I = 8 THEN FILE? = " DRAGN" 

860 IF I = 9 THEN FILE? = "YAWN" 

870 IF I = 10 THEN FILE? = "PITCHMN" 

880 IF I = 11 THEN FILE? = "YAWMN" 

890 IF I = 12 THEN FILE? = "ROLLMN" 

900 REM 
910 REM 

920 GOSUB 1380 'RECORD CALIBRATION DATA 
930 NEXT I 
940 CLS 

950 LOCATE 10, 5: PRINT "CALCULATING CALIBRATION CONSTANTS, 
Figure A. 16 BALCAL . BAS (cont.) 
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PLEASE WAIT" 

960 GOSUB 2100 'CALCULATE K1&K2 FOR ALL PRIME GAGES 
970 GOSUB 2640 'CONVERT COUNTS TO FORCES 
980 GOSUB 3180 'CALCULATE INTERACTION CONSTANTS 
990 REM 

1000 REM STORE CALIBRATION CONSTANTS 
1010 OPEN "C: CONST" FOR OUTPUT AS #1 

1020 WRITE #1 , INCALL, INCALD , INCALY , INCALPM, INCALYM, INCALRM 
1040 WRITE #1 , K1LPOS , K2LPOS , K1DPOS , K2DPOS , K1YPOS , K2YPOS 
1050 WRITE #1 , K1PMPOS , K2PMPOS , K1YMPOS , K2 YMPOS , K1RMPOS , 
K2RMPOS 

1060 WRITE #1 , K1LNEG, K2LNEG,K1DNEG, K2DNEG , K1YNEG, K2YNEG 
1070 WRITE #1,K1PMNEG, K2PMNEG , K1YMNEG , K2YMNEG, K1RMNEG, 
K2RMNEG 

1080 WRITE #1 , DDDL1P, DDDL2P, DYDL1P, DYDL2P, DPMDL1P, DPMDL2P 
, DYMDL1P, DYMDL2P 
1090 WRITE # 1 , DRMDL1 P , DRMDL2 P 

1100 WRITE #1, DLDD1P, DLDD2P, DYDD1P, DYDD2P, DPMDD1P, DPMDD2P 
, DYMDD1P, DYMDD2P 
1110 WRITE #1 , DRMDD1P, DRMDD2P 

1120 WRITE #1, DLDY1P, DLDY2P, DDDY1P, DDDY2P, DPMDY1P, DPMDY2P 
, DYMDY1P, DYMDY2P 
1130 WRITE #1, DRMDY1P, DRMDY2P 

1140 WRITE #1 , DLDPM1P , DLDPM2P , DDDPM1P, DDDPM2P , DYDPM1P , DYD 
PM2P, DYMDPM1P, DYMDPM2P 
1150 WRITE # 1 , DRMDPM1P , DRMDPM2 P 

1160 WRITE #1, DLDYM1P, DLDYM2P, DDDYM1P, DDDYM2P, DYDYM1P, DYD 
YM2P, DPMDYM1P, DPMDYM2P 
1170 WRITE # 1 , DRMDYM1P , DRMDYM2 P 

1180 WRITE #1 , DLDRM1P, DLDRM2P, DDDRM1P, DDDRM2P, DYDRM1P, DYD 
RM2P, DPMDRM1P , DPMDRM2 P 
1190 WRITE #1 , DYMDRM1P, DYMDRM2P 

1200 WRITE #1, DDDL1N, DDDL2N, DYDL1N, DYDL2N, DPMDL1N, DPMDL2N 
, DYMDL1N , DYMDL2N 
1210 WRITE # 1 , DRMDL1N , DRMDL2N 

1220 WRITE # 1 , DLDD1N , DLDD2N , DYDD1N , DYDD2N , DPMDD1N , DPMDD2N 
, DYMDD1N , DYMDD2N 
1230 WRITE # 1 , DRMDD1N , DRMDD2N 

1240 WRITE #1 , DLDY1N, DLDY2N, DDDY1N, DDDY2N, DPMDY1N , DPMDY2N 
, DYMDY1N , DYMDY2N 
1250 WRITE #1 , DRMDY1N, DRMDY2N 

1260 WRITE #1 , DLDPM1N, DLDPM2N, DDDPM1N, DDDPM2N , DYDPM1N , DYD 
PM2N , DYMDPM1N 

1270 WRITE # 1 , DYMDPM2N , DRMDPM1N , DRMDPM2N 

1280 WRITE #1 , DLDYM1N, DLDYM2N, DDDYM1N, DDDYM2N, DYDYM1N, DYD 
YM2N , DPMDYM1N 

1290 WRITE # 1 , DPMDYM2N , DRMDYM1N , DRMDYM2N 

1300 WRITE #1, DLDRM1N, DLDRM2N, DDDRM1N, DDDRM2N, DYDRM1N, DYD 
RM2N , DPMDRM1N 

1310 WRITE # 1 , DPMDRM2N , DYMDRM1N , DYMDRM2N 

Figure A. 16 BALCAL . BAS (cont.) 
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1320 CLOSE #1 

1330 CLS : LOCATE 10,5 

1340 PRINT " CALIBRATION COMPLETE" 

1350 END 
1360 REM 
1370 REM 

13 80 REM RECORD CALIBRATION DATA ■— 

1401 COLOR 15,1: KEY OFF: CLS 

1402 PRINT "CALIBRATION FOR THE LOADING OF THE »;FILE$;" 
COMPONENT" 

1403 PRINT: INPUT "AFTER ALL AMPLIFIERS ARE ZEROED PRESS 
RETURN" ;X 

1404 CLS 

1405 PRINT " ZD ZL ZY ZPM ZYM 

ZRM " 

1406 STAT%=0 : MODE%=0: STOR%=0: COUNT=100: RATE=500 

1407 ADAPT %= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 

1408 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

1409 ZD = 0 : ZPM=0 : ZL=0 : Z YM=0 

1410 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

1411 FOR J = 0 TO 396 STEP 4 

1412 DAT ( J ) = ( DAT% (J)/204.8) -10 

1413 ZD =ZD + DAT ( J ) 

1414 NEXT J 

1415 ZD = ZD/100 

1416 FOR J = 1 TO 397 STEP 4 

1417 DAT ( J) = (DAT% (J)/204.8)-10 

1418 ZL = ZL + DAT ( J) 

1419 NEXT J 

1420 FOR J = 2 TO 398 STEP 4 

1421 DAT ( J) = (DAT% (J)/204.8)-10 

1422 ZPM = ZPM + DAT ( J ) 

1423 NEXT J 

1424 FOR J = 3 TO 399 STEP 4 

1425 DAT ( J) = (DAT% (J)/204.8)-10 

1426 ZYM = ZYM + DAT ( J ) 

1427 NEXT J 

1428 STAT%=0: MODE%=0: STOR%=0 : COUNT=100: RATE=500 

1429 ADAPT%= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=1 

1430 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR%, COUNT, RATE, DAT1% (0) ,STAT%) 

1431 ZY =0 : ZRM=0 

1432 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

1435 FOR J = 0 TO 198 STEP 2 

1436 DAT1 ( J) = (DAT1% (J)/204.8)-10 

1437 ZRM = ZRM + DAT1(J) 

1438 NEXT J 



Figure A. 16 BALCAL . BAS (cont.) 
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1439 ZRM = ZRM/100 

1440 FOR J = 1 TO 199 STEP 2 

1441 DAT1 ( J) = (DAT1% (J)/204.8)-10 

1442 ZY = ZY + DAT1(J) 

1443 NEXT J 

1444 ZL=ZL/100 : ZPM=ZPM/100 : ZYM=ZYM/100 : ZY=ZY/100 

1445 LOCATE 3,1: PRINT USING "+#.###";ZD: LOCATE 3,10: PRI 
NT USING "+#.###" ;ZL 

1446 LOCATE 3,19: PRINT USING "+#.###" ;ZY 

1447 LOCATE 3,28: PRINT USING "+#.###" ;ZPM 

1448 LOCATE 3,37: PRINT USING "+#.###»; ZYM 

1449 LOCATE 3,46: PRINT USING "+#.###"; ZRM 

1451 REM 

1452 PRINT -.PRINT 

1453 INPUT "AFTER PLACING ALL CAL SWITCHES TO + SETTING HIT 
RETURN" ;X 

1454 PRINT 

1455 PRINT " CALD CLL CALY CALPM CALYM 

CALRM" 

1456 STAT%=0: MODE%=0 : STOR%=0: COUNT=100: RATE=500 

1457 ADAPT %= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 

1458 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% ,MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

1459 CALD=0 : CLL=0 : CALYM=0 : CALPM=0 

1460 IF STAT%<> 0 , THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

1461 FOR J = 0 TO 396 STEP 4 

1462 DAT ( J) = ( DAT% (J)/204.8)-10 

1463 CALD =CALD + DAT ( J) 

1464 NEXT J 

1465 CALD = CALD/ 100 

1466 FOR J = 1 TO 397 STEP 4 

1467 DAT ( J) = (DAT% (J)/204.8)-10 

1468 CLL = CLL + DAT ( J) 

1469 NEXT J 

1470 FOR J = 2 TO 398 STEP 4 

1471 DAT ( J) = ( DAT% (J)/204.8)-10 

1472 CALPM = CALPM + DAT ( J) 

1473 NEXT J 

1474 FOR J = 3 TO 399 STEP 4 

1475 DAT ( J) = (DAT% (J)/204.8)-10 

1476 CALYM = CALYM + DAT ( J ) 

1477 NEXT J 

1478 STAT%=0 : MODE%=0: STOR%=0: COUNT=100: RATE=500 

1479 ADAPT%= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=1 

1480 CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

1481 CALY =0 : CALRM=0 

1482 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###»; 
STAT% : END 
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1485 FOR J = 0 TO 198 STEP 2 

1486 DAT1 ( J) = (DAT1% (J)/204.8)-10 

1487 CALRM = CALRM + DAT1(J) 

1488 NEXT J 

1489 CALRM = CALRM/100 

1490 FOR J = 1 TO 199 STEP 2 

1491 DATl(J)=(DATl%(J)/204.8)-10 

1492 CALY = CALY + DAT1 ( J) 

1493 NEXT J 

1494 CLL=CLL/100 : CALPM=CALPM/ 100 : CALYM=CALYM/ 100 : CALY=CALY 
/100 

1495 LOCATE 10,1: PRINT USING "+#.###"; CALD 

1496 LOCATE 10,10: PRINT USING "+#.###" ;CLL 

1497 LOCATE 10,19: PRINT USING "+#.###" ;CALY 

1498 LOCATE 10,28:. PRINT USING "+#.###"; CALPM 

1499 LOCATE 10,37: PRINT USING "+#.###" ;CALYM 

1500 LOCATE 10,46: PRINT USING "+#.###"; CALRM 

1501 REM 

1502 REM 



1503 


IF 


I 


= 


l 


THEN 


INCALL = 


CLL - : 


SL 


1510 


IF 


I 


= 


2 


THEN 


INCALD = 


CALD - 


ZD 


1520 


IF 


I 


=r 


3 


THEN 


INCALY = 


CALT - 


ZY 


1530 


IF 


I 


= 


4 


THEN 


INCALPM 


= CALPM 


- ZPM 


1540 


IF 


I 


= 


5 


THEN 


INCALYM 


= CALYM 


- ZYM 


1550 


IF 


I 


= 


6 


THEN 


INCALRM 


= CALRM 


- ZRM 



1611 PRINT: INPUT "PRESS ENTER (RETURN) TO CONTINUE" ;X 
1620 CLS: PRINT "RETURN ALL CAL SWITCHES TO CENTER POSITION" 
1630 PRINT: 

1640 PRINT "LOAD THE ";FILE$;" GAGE FROM 0 TO 20 POUNDS AND 
THEN BACK TO 0" 

1650 PRINT "POUNDS IN 1 POUND INCREMENTS. THERE SHOULD BE 
TWO READINGS FOR" 

1660 PRINT "EACH POUND WEIGHT EXCEPT FOR THE LOAD AT 20 
POUNDS WHICH WILL HAVE" 

1670 PRINT "ONLY ONE READING" 

1680 PRINT 

1690 PRINT "PRESS THE FI KEY WHEN THE LOADING IS FINISHED" 
1700 PRINT "PRESS THE F2 KEY WHEN READY TO RECORD THE DATA 
FOR THAT LOAD" 



1701 


PRINT: 


INPUT 


"PRESS ENTER 


(RETURN) 


TO CONTINUE" ;X 


1702 


CLS 












1703 


PRINT 


" DRAG 


LIFT 


YAW 


PITCH 


YAW 




ROLL 


ii 










1704 


PRINT 


ii 






MOM. 


MOM. 



MOM. " 

1710 SOAP = 0: N=4 
1720 FOR K = 1 TO 100 

1730 ON KEY ( 1) GOSUB 1830 'SET STOP FLAG 

1740 ON KEY ( 2 ) GOSUB 1841 'RECORD DATA 

1750 KEY ( 1) ON: KEY (2) ON 
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1760 IF SOAP = 2 THEN GOTO 1780 

1770 IF SOAP = 1 THEN GOTO 1920 

1775 GOTO 1730 

1780 SOAP = 0 

1790 NOD% = K 

1800 NEXT K 

1810 GOTO 1920 

1830 REM SET STOP FLAG 

1838 SOAP = 1 

1839 RETURN 

1840 REM STEPS TO RECORD DATA ' 

1841 STAT%=0 : MODE%=0: STOR%=0: COUNT=100: RATE=500 

1842 ADAPT %= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 

1843 CALL AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT%) 

1844 D (K) =0 : L(K) =0 : YM (K) =0 : PM (K) =0 

1845 IF STAT%<> 0 , THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

1846 FOR J = 0 TO 396 STEP 4 

1847 DAT ( J) = (DAT% (J)/204.8)-10 

1848 D (K) =D(K) + DAT ( J ) 

1849 NEXT J 

1850 D (K) = D(K)/100 

1851 FOR J = 1 TO 397 STEP 4 

1852 DAT (J) = (DAT% (J)/204.8)-10 

1853 L(K) = L(K) + DAT ( J) ' 

1854 NEXT J 

1855 FOR J = 2 TO 398 STEP 4 

1856 DAT (J) = (DAT% (J)/204.8)-10 

1857 PM (K) = PM (K) + DAT ( J ) 

1858 NEXT J 

1859 FOR J = 3 TO 399 STEP 4 

1860 DAT ( J ) = ( DAT% (J)/204.8)-10 

1861 YM(K) = YM (K) + DAT ( J ) 

1862 NEXT J 

1863 STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 

1864 ADAPT %= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=1 

1865 CALL AINSC ( ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT ; RATE , DAT1% ( 0 ) , STAT% ) 

1866 Y (K) =0 : RM (K) =0 

1867 IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

1868 FOR J = 0 TO 198 STEP 2 

1869 DAT1 ( J) = ( DAT1% (J)/204.8)-10 

1870 RM (K) = RM (K) + DAT1(J) 

1871 NEXT J 

1872 RM (K) = RM(K)/100 

1873 FOR J = 1 TO 199 STEP 2 

1874 DAT1 (J) = (DAT1% (J)/204.8)-10 

1875 Y (K) = Y(K) + DAT1(J) 
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1876 

1877 

1878 
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1933 

1934 

1935 
1940 
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1970 

1971 
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1980 

1981 

1982 

1983 

1984 
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NEXT J 

L(K)=L(K)/100:PM(K)=PM(K)/100:YM(K) =YM (K) /100 : Y (K) =Y ( 
K)/100 

IF N> 23 THEN N=4 : CLS : PRINT " DRAG LIFT YAW 

PITCH 
MOM. 

LOCATE 

LOCATE N,10: PRINT 
LOCATE N, 19 : PRINT 
LOCATE N, 28 : PRINT 
LOCATE N,37: PRINT 
LOCATE N, 46 : PRINT 
N=N+1 



23 THEN N=4 
YAW 
MOM. 
N, 1: PRINT 



CLS: PRINT " DRAG 
ROLL "-.PRINT " 

MOM. " 

USING "+#.###" ;D(K) 
USING 
USING 
USING 
USING 
USING 



"+# . ###" ;L(K) 
"+#.###" ;Y(K) 
"+#.###"; PM (K) 
»+#.###" ;YM(K) 
"+#.###";RM(K) 



SOAP = 2 

RETURN 'GO BACK FOR NEXT DATA POINT 
i CORRECT FOR DRIFT 

DIFFL =(L(NOD%) -L(l) ) : DIFFD= (D (NOD%) -D ( 1) ) : DIFFY= ( Y ( 
NOD%) -Y ( 1) ) 

DIFFPM= ( PM (NOD%) -PM ( 1) ) : DIFFYM= ( YM (NOD% ) -YM ( 1) ) : DIFFR 
M=(RM(NOD%) -RM( 1) ) 

D=N0D%-1 



CORL=DIFFL/D : CORD=DIFFD/D : CORY=DIFFY/D : CORPM=DIFFPM/D 

: CORYM=DIFFYM/D : CORRM=DIFFRM/D 

A=1 



FOR K = 2 TO NOD% 

L (K) =L (K) — (A*CORL) 

D ( K) =D (K) - (A*CORD) 

Y(K)=Y(K)-(A*CORY) 

PM (K) =PM (K) - (A*CORPM) 

YM ( K) =YM ( K) - ( A*CORYM) 

RM(K) =RM (K) - (A*CORRM) 

A = A+l 
NEXT K 

' WRITE DATA TO FILE 

FILE$ = "C:"+FILE$ 

OPEN FILE$ FOR OUTPUT AS #1 
WRITE #1, ZL, ZD, ZY, ZPM, ZYM, ZRM, NOD% 

WRITE # 1 , CLL , CALD , CALY , CALPM , CALYM , CALRM 
IF 1=1 OR 1=2 OR 1=3 THEN GOTO 2079 
IF 1=7 OR 1=8 OR 1=9 THEN GOTO 2067 
LBS = 0 

FOR J = 1 TO 11 

IF 1=4 OR 1=10 THEN Dl=20 

IF 1=5 OR 1=11 THEN Dl=4 . 5 

IF 1=6 OR 1=12 THEN Dl=11.5 

MOM = LBS*D1 'CONVERT TO MOMENTS 

IF 1=10 OR 1=11 OR 1=12 THEN MOM=MOM*(-l) 

WRITE #1 , L( J) , D ( J) , Y ( J) , PM ( J ) ,YM(J) ,RM(J) ,MOM 

LBS = LBS+1 

NEXT J 
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1986 LBS = 9 

1987 FOR J = 12 TO NOD% 

1988 IF 1=4 OR 1=10 THEN Dl=20 

1989 IF 1=5 OR 1=11 THEN Dl=4 . 5 

1990 IF 1=6 OR 1=12 THEN Dl-11.5 

1991 MOM = LBS*D1 'CONVERT TO MOMENTS 

1992 IF 1=10 OR 1=11 OR 1=12 THEN MOM=MOM*(-l) 

1993 WRITE #1,L(J) ,D(J) ,Y(J) ,FM(J) ,YM(J) ,RM(J) ,MOM 

1994 LBS = LBS-1 

1995 NEXT J 

1996 CLOSE #1 

1997 RETURN 'GO BACK TO SET UP NEXT CALIBRATION 

2067 LBS = 0 

2068 FOR J = 1 TO 21 

2069 WRITE #1 , L ( J) , D ( J) , Y ( J) , PM ( J) , YM ( J) , RM ( J) , LBS 

2070 LBS = LBS - 1 

2071 NEXT J 

2072 LBS = -19 

2073 FOR J = 22 TO NOD% 

2074 WRITE # 1 , L ( J) , D ( J) , Y ( J) , PM ( J) , YM ( J) , RM ( J) , LBS 

2075 LBS = LBS + 1 

2076 NEXT J 

2077 CLOSE #1 

2078 RETURN 'GO BACK TO SET UP NEXT CALIBRATION 

2079 LBS = 0 

2080 FOR J = 1 TO 21 

2081 WRITE #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) , LBS 

2082 LBS = LBS + 1 

2083 NEXT J 

2084 LBS = 19 

2085 FOR J = 22 TO NOD% 

2086 WRITE # 1 , L ( J) , D ( J) , Y ( J) , PM ( J) , YM ( J) , RM ( J) , LBS 

2087 LBS = LBS - 1 

2088 NEXT J 

2089 CLOSE #1 

2090 RETURN 'GO BACK TO SET UP NEXT CALIBRATION 

2100 REM 

2110 REM CALCULATE K1&K2 FOR PRIME GAGES 



2120 


FOR I 


— 


1 


TO 12 




2130 


IF 


I 


= 


1 


THEN 


FILE$ = "LIFTP'' 


2140 


IF 


I 


= 


2 


THEN 


FILE$ = "DRAGP" 


2150 


IF 


I 


= 


3 


THEN 


FILE$ = "YAWP" 


2160 


IF 


I 


= 


4 


THEN 


FILE$ = "PITCHMP" 


2170 


IF 


I 


= 


5 


THEN 


FILE$ = " YAWMP" 


2180 


IF 


I 


= 


6 


THEN 


FILE$ = "ROLLMP" 


2190 


IF 


I 


= 


7 


THEN 


FILE$ = "LIFTN" 


2200 


IF 


I 


= 


8 


THEN 


FILE$ = "DRAGN" 


2210 


IF 


I 


= 


9 


THEN 


FILE$ = "YAWN" 


2220 


IF 


I 


= 


10 THEN FILE$ = "PITCHMN 


2230 


IF 


I 


= 


11 THEN FILE$ = " YAWMN" 
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IF I = 12 THEN FILE$ = "ROLLMN" 

GOSUB 2295 'READ FILES AND PERFORM CALCULATIONS 
NEXT I 

RETURN 'GO BACK TO CONVERT COUNTS 

REM 

REM 

A=0 : B=0 : C=0 : D=0 : E=0 : 

OPEN "C:"+FILE$ FOR INPUT AS #2 
INPUT #2 , ZL, ZD, ZY , ZPM, ZYM, ZRM, NOD% 

INPUT # 2 , CLL , CALD , CALY , CALPM , CALYM , CALRM 
FOR J = 1 TO NOD% 

INPUT #2 , L ( J) , D ( J) , Y ( J) , PM ( J ) ,YM(J) ,RM(J) ,LBS(J) 



IF 


I 


= 


1 


OR 


I 


= 


7 THEN X = 


L(J) 


IF 


I 


= 


2 


OR 


I 


= 


8 THEN X = 


D ( J) 


IF 


I 


= 


3 


OR 


I 


= 


9 THEN X = 


Y(J) 


IF 


I 


= 


4 


OR 


I 


= 


10 THEN X 


= PM ( J) 


IF 


I 


= 


5 


OR 


I 


= 


11 THEN X 


= YM ( J) 


IF 


I 


= 


6 


OR 


I 


= 


12 THEN X 


= RM ( J ) 



A = A + (X A 2 ) 

B = B + (X A 3 ) 

C = C + (X A 4 ) 

D = D + (X*LBS ( J) ) 

E = E + (X*X*LBS ( J) ) 

NEXT J 
CLOSE #2 

K2 = ( (D/A) - (E/B) ) / ( (B/A) - (C/B) ) 
K1 = (D/A) - (K2* (B/A) ) 



IF 


I 


= 


1 


THEN 


K1LPOS = 


Kl: 


K2LPOS = 


K2 


IF 


I 


= 


2 


THEN 


K1DPOS = 


Kl: 


K2DPOS = 


K2 


IF 


I 


= 


3 


THEN 


K1YPOS = 


Kl: 


K2YPOS = 


K2 


IF 


I 


= 


4 


THEN 


K1PMPOS : 


= Kl 


: K2PMPOS 


= K2 


IF 


I 


= 


5 


THEN 


K1YMPOS = 


= Kl 


: K2YMPOS 


= K2 


IF 


I 


= 


6 


THEN 


K1RMPOS = 


= Kl 


: K2RMPOS 


= K2 


IF 


I 


= 


7 


THEN 


K1LNEG = 


Kl: 


K2LNEG = 


K2 


IF 


I 


= 


8 


THEN 


K1DNEG = 


Kl: 


K2DNEG = 


K2 


IF 


I 


= 


9 


THEN 


K1YNEG = 


Kl: 


K2YNEG = 


K2 



IF I = 10 THEN K1PMNEG = K1 : K2PMNEG = K2 

IF I = 11 THEN K1YMNEG = Kl: K2YMNEG = K2 

IF I = 12 THEN K1RMNEG = Kl: K2RMNEG = K2 

RETURN 'GO BACK TO CALCULATE K1&K2 FOR NEXT FILE 

REM 

REM CONVERT COUNTS TO FORCES 



FOR I 


= 


1 


TO 12 








IF 


I 


= 


1 


THEN 


FILE$ 


= 


"LIFTP" 


IF 


I 


= 


2 


THEN 


FILE$ 


= 


" DRAGP" 


IF 


I 


= 


3 


THEN 


FILE$ 


= 


"YAWP" 


IF 


I 


= 


4 


THEN 


FILE$ 


= 


"PITCHMP" 


IF 


I 


= 


5 


THEN 


FILE$ 


= 


"YAWMP" 


IF 


I 


= 


6 


THEN 


FILE$ 


= 


"ROLLMP" 


IF 


I 


= 


7 


THEN 


FILE$ 


= 


" LIFTN" 
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2740 


IF 


I 


= 


8 


THEN FILE$ = 


" DRAGN" 


2750 


IF 


I 


= 


9 


THEN FILE$ = 


"YAWN" 


2760 


IF 


I 


= 


10 


THEN FILES 


= "PITCHMN" 


2770 


IF 


I 


= 


11 


THEN FILES 


= " YAWMN" 


2780 


IF 


I 


= 


12 


THEN FILES 


= "ROLLMN" 



2790 GOSUB 2835 'READ IN FILE AND CONVERT 
2800 NEXT I 

2810 RETURN 'GO BACK TO CALCULATE INTERACTON CONSTANTS 
2820 REM 
2830 REM 

2835 A=0:B=0:C=0:D=0:E=0:F=0 

2840 OPEN "C:"+FILE$ FOR INPUT AS #2 

2850 INPUT #2 , ZL, ZD , ZY , ZPM, ZYM, ZRM, NOD% 

2 860 INPUT # 2 , CLL , CALD , CALY , CALPM , CALYM , CALRM 
2870 FOR J = 1 TO NOD% 

2880 INPUT #2,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 
2890 A=( (INCALL/CLL-ZL) * (L(J) -ZL) ) 

2900 B= ( (INCALD/CALD-ZD) * (D ( J) -ZD) ) 

2910 C=( (INCALY/CALY-ZY) * (Y(J) -ZY) ) 

2920 D=( (INCALPM/ CALPM- ZPM) *(PM(J) -ZPM) ) 

2930 E=( (INCALYM/ CALYM- ZYM) *(YM(J) -ZYM) ) 

2940 F= ( (INCALRM/ CALRM-ZRM) * (RM ( J) -ZRM) ) 

2950 IF L ( J) < 0 THEN K1 = K1LNEG: K2 = K2LNEG ELSE K1 
= K1LPOS : K2 = K2LPOS 
2960 L(J) = (K1*A) + (K2 * (A A 2 ) ) 

2970 IF D ( J) <0 THEN K1 = K1DNEG: K2 = K2DNEG ELSE K1 
= K1DPOS : K2 = K2DPOS 
2980 D ( J) = (K1*B) + (K2 * ( B A 2 ) ) 

2990 IF Y ( J) <0 THEN K1 = K1YNEG: K2 = K2YNEG ELSE K1 
= K1YPOS : K2 = K2YPOS 
3000 Y(J)=(K1*C) + (K2* (C A 2 ) ) 

3010 IF PM ( J) < 0 THEN K1=K1PMNEG : K2=K2 PMNEG ELSE K1=K1 
PMPOS : K2=K2 PMPOS 

3020 PM( J) = (K1*D) + (K2* (D A 2) ) 

3030 IF YM ( J) < 0 THEN K1=K1YMNEG : K2=K2 YMNEG ELSE K1=K1 
YMPOS : K2=K2YMPOS 

3040 YM ( J) = ( K1*E) + (K2*(E A 2)) 

3050 IF RM ( J) < 0 THEN K1=K1RMNEG : K2=K2RMNEG ELSE K1=K1 
RMPOS : K2=K2RMPOS 

3060 RM ( J) = (K1*F) + (K2*(F A 2)) 

3070 NEXT J 

3080 CLOSE #2 

3090 REM SAVE CONVERTED COUNTS 

3100 OPEN "B : "+FILE$ FOR OUTPUT AS #1 

3110 WRITE # 1 , ZL, ZD , ZY , ZPM, ZYM , ZRM , NOD% 

3120 WRITE #1, CLL, CALD, CALY, CALPM, CALYM, CALRM 
3130 FOR J = 1 TO NOD% 

3140 WRITE #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 

3150 NEXT J 

3155 CLOSE #1 
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3160 RETURN 'GO BACK TO CALCULATE INTERACTION CONSTANTS 

3170 REM 

3180 REM CALCULATE INTERACTION CONSTANTS 
3190 FOR I = 1 TO 12 

3200 IF I = 1 THEN FILE? = " LI FTP" 

3210 IF I = 2 THEN FILE? = "DRAGP" 

3220 IF I = 3 THEN FILE? = "YAWP" 

3230 IF I = 4 THEN FILE? = "PITCHMP" 

3240 IF 1=5 THEN FILE? = "YAWMP" 

3250 IF I = 6 THEN FILE? = "ROLLMP" 

3260 IF I = 7 THEN FILE? = "LIFTN" 

3270 IF 1=8 THEN FILE? = "DRAGN" 

3280 IF 1=9 THEN FILE? = "YAWN" 

3290 IF I = 10 THEN FILE? = "PITCHMN" 

3300 IF I = 11 THEN FILE? = "YAWMN" 

3310 IF I = 12 THEN FILE? = "ROLLMN" 

3320 GOSUB 3365 'READ IN FILES AND CALCULATE CONSTANTS 

3330 NEXT I 

3340 RETURN 'GO BACK TO SAVE CALIBRATION DATA 

3350 REM 
3360 REM 

3365 A=0 : B=0 : C=0 : D1=0 : D2=0 : D3=0 : D4=0 : D5=0 : E1=0 : E2=0 : E3=0 : E 
4=0 : E5=0 

3370 OPEN "B: "+FILE? FOR INPUT AS #2 
3380 INPUT #2,ZL,ZD,ZY, ZPM, ZYM, ZRM,NOD% 

3 39 0 INPUT # 2 , CLL , CALD , CALY , CALPM , CALYM , CALRM 
3410 FOR J = 1 TO NOD% 

3420 INPUT #2,L(J) ,D(J) , Y(J) ,PM(J) , YM(J) ,RM(J) ,LBS(J) 
3430 IF 1=1 OR 1=7 THEN X=L(J) :Y1=D(J) :Y2=Y(J) :Y3=PM(J) 

: Y4=YM ( J) :Y5=RM(J) 

3440 IF 1=2 OR 1=8 THEN X=D ( J) : Y1=L ( J) : Y2=Y ( J) : Y3=PM ( J) 

: Y4=YM ( J) : Y5=RM ( J) 

3450 IF 1=3 OR 1=9 THEN X=Y ( J) : Y1=L ( J) : Y2=D ( J) : Y3=PM ( J) 

: Y4=YM(J) : Y5=RM(J) 

3460 IF 1=4 OR 1=10 THEN X=PM ( J) : Y1=L ( J) : Y2=D ( J) : Y3=Y ( J) 
: Y4=YM( J) : Y5=RM ( J) 

3470 IF 1=5 OR 1=11 THEN X=YM ( J) : Y1=L( J) : Y2=D ( J) : Y3=Y ( J) 
:Y4=PM(J) : Y5=RM(J) 

3480 IF 1=6 OR 1=12 THEN X=RM ( J) : Y1=L ( J) : Y2=D ( J) : Y3=Y ( J) 
: Y4=PM(J) :Y5=YM(J) 



3490 


A = 


< A + 


(X 


A 2 ) 


3500 


B = 


> B + 


(X 


-3) 


3510 


C = 


> c + 


(X 


A 4 ) 


3520 


D1 


= D1 


+ 


(X*Y1) 


3530 


El 


= El 


+ 


( (X A 2 ) *Y1) 


3540 


D2 


= D2 


+ 


(X*Y2 ) 


3550 


E2 


= E2 


+ 


( (X A 2 ) *Y2 ) 


3560 


D3 


= D3 


+ 


(X*Y3 ) 


3570 


E3 


= E3 


+ 


( (X A 2) *Y3 ) 


3580 


D4 


= D4 


+ 


(X*Y4 ) 
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E4 = E4 + ( (X A 2 ) *Y4 ) 

D5 = D5 + (X*Y5 ) 

E5 = E5 + ( (X A 2 ) *Y5) 

NEXT J 
CLOSE #2 

K12 = ( (D1/A)-(E1/B) )/((B/A)-(C/B) ) 

Kll = (Dl/A) - (K12* (B/A) ) 

K22 = ( (D2/A) - (E2/B) ) / ( (B/A) - (C/B) ) 

K21 = (D2/A) - (K22* (B/A) ) 

K32 = ( (D3/A)-(E3/B) ) / ( (B/A) — (C/B) ) 

K31 = (D3/A) - (K32* (B/A) ) 

K42 = ( (D4/A) -(E4/B) )/ ( (B/A) - (C/B) ) 

K41 = (D4/A) - (K42* (B/A) ) 

K52 = ( (D5/A) -(E5/B) )/ ( (B/A) - (C/B) ) 

K51 = (D5/A) - (K52* (B/A) ) 

IF 1=1 THEN DDDL1P=K11:DDDL2P=K12:DYDL1P=K21:DYDL2P=K 
22 : DPMDL1P=K3 1 

IF 1=1 THEN DPMDL2 P=K3 2 : DYMDL1P=K4 1 : D YMDL2 P=K4 2 : DRMDL 
1P=K5 1 : DRMDL2 P=K5 2 

IF 1=2 THEN DLDD1P=K11:DLDD2P=K12:DYDD1P=K21:DYDD2P=K 
22 :DPMDD1P=K31 

IF 1=2 THEN DPMDD2P=K3 2 : DYMDD1P=K41 : DYMDD2P=K42 : DRMDD 
1P=K51 : DRMDD2P=K52 

IF 1=3 THEN DLDY1P=K11:DLDY2P=K12:DDDY1P=K21:DDDY2P=K 
22 : DPMDY1P=K3 1 

IF 1=3 THEN DPMDY2P=K32 : DYMDY1P=K41 : DYMDY2P=K42 : DRMDY 
1P=K5 1 : DRMDY2 P=K5 2 

IF 1=4 THEN DLDPM1P=K11:DLDPM2P=K12 :DDDPM1P=K21:DDDPM 
2P=K22 : DYDPM1P=K3 1 

IF 1=4 THEN DYDPM2P=K32 : DYMDPM1P=K41 : DYMDPM2P=K42 : DRM 
DPM1P=K51 : DRMDPM2P=K52 

IF 1=5 THEN DLDYM1P=K11:DLDYM2P=K12 : DDDYM1P=K2 1 : DDDYM 
2P=K2 2 : DYDYM1P=K3 1 

IF 1=5 THEN DYDYM2 P=K3 2 : DPMDYM1P=K4 1 : DPMDYM2 P=K4 2 : DRM 
DYM1P=K51 : DRMD YM2 P=K5 2 

IF 1=6 THEN DLDRM1P=K11 : DLDRM2P=K12 : DDDRM1P=K2 1 : DDDRM 
2P=K22 : DYDRM1P=K3 1 

IF 1=6 THEN DYDRM2 P=K3 2 : DPMDRM1P=K4 1 : DPMDRM2 P=K4 2 : DYM 
DRM1P=K5 1 : DYMDRM2 P=K5 2 

IF 1=7 THEN DDDL1N=K11:DDDL2N=K12:DYDL1N=K21:DYDL2N=K 
22 : DPMDL1N=K3 1 

IF 1=7 THEN DPMDL2N=K32:DYMDL1N=K41:DYMDL2N=K42: DRMDL 
1N=K51 : DRMDL2N=K52 

IF 1=8 THEN DLDD1N=K11:DLDD2N=K12:DYDD1N=K21:DYDD2N=K 
22 : DPMDD1N=K3 1 

IF 1=8 THEN DPMDD2N=K3 2 : DYMDD1N=K4 1 : DYMDD2N=K4 2 : DRMDD 
1N=K51 : DRMDD2N=K52 

IF 1=9 THEN DLDY1N=K11:DLDY2N=K12:DDDY1N=K21:DDDY2N=K 
22 : DPMDY1N=K3 1 

IF 1=9 THEN DPMDY2N=K3 2.*DYMDY1N=K41:DYMDY2N=K42: DRMDY 
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3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 



1N=K51 : DRMDY2N=K52 

IF 1=10 THEN DLDPM1N=K11:DLDPM2N=K12:DDDPM1N=K21:DDDP 
M2N=K22 : DYDPM1N=K31 

IF 1=10 THEN DYDPM2N=K32:DYMDPM1N=K41:DYMDPM2N=K42:DR 
MDPM1N=K51 

IF 1=10 THEN DRMDPM2N=K52 

IF 1=11 THEN DLDYM1N=K11 : DLDYM2N=K12 : DDDYM1N=K2 1 : DDDY 
M2N=K22 : DYDYM1N=K31 

IF 1=11 THEN DYDYM2N=K32:DPMDYM1N=K41:DPMDYM2N=K42:DR 
MDYM1N=K51 

IF 1=11 THEN DRMDYM2N=K52 

IF 1=12 THEN DLDRM1N=K11:DLDRM2N=K12:DDDRM1N=K21:DDDR 
M2N=K22 : DYDRM1N=K31 

IF 1=12 THEN DYDRM2N=K32 :DPMDRM1N=K41:DPMDRM2N=K42 :DY 
MDRM1N=K51 

IF 1=12 THEN D YMDRM2 N=K5 2 

RETURN 'GOBACK TO SAVE CALIBRATION DATA 
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APPENDIX B 



TABLES 



INITIAL 


SET-UP FOR 


THE MODEL 


8255 TRANSDUCER 


AMPLIFIER 


Amp . # 


Component 


Gain Set. 


Filter Set. 


EX. Volt. 


1 


DRAG 


MAX VAR. 


1 


+ 5.0 


2 


LIFT 


MAX VAR. 


1 


+ 6 . 5 


3 


PITCH M. 


IK 


1 


+ 5.0 


4 


YAW M. 


MAX VAR. 


1 


+ 5.0 


5 


ROLL M. 


MAX VAR. 


1 


+ 5.0 


6 


YAW 


MAX VAR. 


1 


+ 5.0 


7 


AOA 


IK 


1 


+ 5.0 
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Table B.13 Loading for Negative Rolling Moment 
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